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We study the shape, elasticity and fluctuations of the recently predictedS and subsequently ob- 
served (in numerical simulations)!! tubule phase of anisotropic membranes, as well as the phase 
transitions into and out of it. This novel phase lies between the previously predicted flat and crum- 
pled phases, both in temperature and in its physical properties: it is crumpled in one direction, 
and extended in the other. Its shape and elastic properties are characterized by a radius of gyra- 
tion exponent v and an anisotropy exponent z. We derive scaling laws for the radius of gyration 
Rg{L±, L y ) (i.e. the average thickness) of the tubule about a spontaneously selected straight axis 
and for the tubule undulations h rms (L±, L y ) transverse to its average extension. We show that for 
square membranes (with intrinsic size L± — L y = Li), Rg oc L", and h rms oc L 1 ~ v " z ^ 2 , w ith r/ K 
a bending rigidity anomalous elasticity exponent related to v and z. For phantom (i.e. non-self- 
avoiding) membranes, we predict v = 1/4, z — 1/2 and n K = 0, exactly, in excellent agreement with 
simulations. For D = 2 dimensional membranes embedded in the space of dimension d < 11, self- 
avoidance greatly swells the tubule and suppresses its wild transverse undulations, changing its shape 
exponents u, z, and rj K . For a D-dimensional membrane embedded in d > d* (d*(D = 2) > 7/2), 
n K = and z = (£>-l+2v)/3, while for d < d*, rj K > and z = (Z?-l+2i/)/(3-r/«). "Flory" theory 
yields, in the physical case of D = 2 and d = 3, v = 3/4, while the recent 11 — e expansion resultsa 
yields v = 0.52. The actual value of v probably lies closer to the Flory estimate, between these two 
limits. We give detailed scaling results for the shape of the tubule of an arbitrary aspect ratio, i.e. 
for the tubule thickness, its transverse undulations, and a variety of other correlation functions, as 
well as for the anomalous elasticity of the tubules, in terms of v and z. Finally we present a scaling 
theory for the shape and specific heat near the continuous transitions into and out of the tubule 
phase and perform detailed renormalization group calculations for the crumpled-to-tubule transition 
for phantom membranes. 
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I. INTRODUCTION 



Tethered membraneaj became a subject of great in- 
terest when it was theoretically predictedQ that, unlike 
polymers, which are always orientationally disordered, 
membranes can exhibit two distinct phases: crumpled 
and flat, with a "crumpling" transition between them. 
The fiat phase is particularly novel and intriguing, be- 
cause it provides an example of a two dimensional system 
with a continuous symmetry that nonetheless exhibits a 
long-ranged order (specifically, long-ranged orientational 
order in the normal to the membrane) in appareat viola- 
tion of the Hohenberg-Mermin- Wagner theoremE 
ordering is made possible by "anomalous elasticity"!! 
thermal fluctuations infinitely enhance the bending rigid- 
ity k of the membrane at long wavelengths, thereby sta- 
bilizing the orientational order against these very fluc- 
tuations. This is perhaps the most dramatic illustration 
yet found of the phenomenon of "order from disorder" . 

Rich as these phenomena are, most past theoretical 



worku has been restricted to isotropic membranes. In a 
recent paperEI we extended these considerations to intrin- 
sically anisotropic membranes (e.g., polymerized mem- 
branes with in-plane tilt orderLI) and found, astonishingly, 
that anisotropy, a seemingly innocuous generalization, 
actually leads to a wealth of new phenomena. Most dra- 
matically, we found an entire new phase of membranes, 
which we called the "tubule" phase, ubiquitously inter- 
venes between the high temperature crumpled and low 
temperature "flat" phases. The defining property of the 
tubule phase is that it is crumpled in one of the two 
membrane directions, but "flat" (i.e., extended) in the 
other. Its average shape is a long, thin cylinder of length 
R y = L y x 0(1) and radius Ra{L±) <C L±, where L y 
and L± are the dimensions the membrane would have in 
the extended and crumpled directions respectively, were 
it to be flattened out. It should be clarified here that we 
use the term "cylinder" extremely loosely; as illustrated 
in Fig.^, a cross section of the membrane perpendicular 
to the tubule axis (y) will look as disordered as a flexible 
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polymer. These tubules, occurring as a low tempera- 
ture phase of anisotropic polymerized membranes, have 
little in common (and therefore should not be confused) 
with micro-rtubules that are found in liquid phospholipid 
membraneglJ. 

Only in—the special case of perfectly isotropic 
membraneslid is it possible for the membrane to undergo a 
direct transition from the flat to the crumpled phase. The 
theoretically predicted!! and recently observedo phase di- 
agram is shown in Fig|l]. 



1- TUBULE 
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FIG. 1. Phase diagram for anisotropic tethered membranes 
showing the new tubule and previously studied flat and crum- 
pled phases. 
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with such order have already been founcEj'Eij; it should be 
possible to polymerize these without destroying the tilt 
order. Secondly, membranes could be fabricated by cross .- 
linking DNA molecules trapped in a fluid membranec3E3. 
Performing the cross-linking in an applied electric field 
would align the DNA and "freeze in" the anisotropy in- 
duced by the electric field, which could then be removed. 

The tubule cross-sectional radius Rq, ( hereafter called 
the radius of gyration) , and its undulations h rma trans- 
verse to its average axis of orientation, obey the scaling 
laws: 



Rg(L±, Ly) 

hrms(L±, Ly) 

where Q = v/z, 



= L y S h (L y /Ll) 



1 



-(1 + 2!/) , 



(1.1) 
(1.2) 



(1.3) 



we have specialized in Eq.L3 to D = 2 (with general 
expression for a D-dimensional membrane given in the 
main text), the universal exponents v and z are < 1, r] K 
is the anomalous elasticity exponent for the tubule bend- 
ing rigidity n (as defined by n ~ L^ K , also see below), and 
for convenience we chose to measure the intrinsic lengths 
L y and L±_ in units of the ultraviolet cutoff, set approxi- 
mately by the the monomer (e.g. phospholipid) size. 



The direct crumpling transition studied previously oc- 
curs in our more generic model only for that special set 
of cuts through the phase diagram (like P2) that pass 
through the origin. Generic paths (like Pi) will ex- 
perience two phase transitions, crumpled-to-tubule, and 
tubule-to-flat, that are in new, heretofore uninvestigated 
universality classes. 

This prediction was recently dramatically confirmed 
in Monte Carlo simulations of phantom (i.e., non-self- 
avoiding) membranes by Bowick, Falcioni and Thorleif- 
sson (BFT).cl. They simulated membranes with differ- 
ent bare bending moduli n x and n y in the orthogonal x 
and y directions. As temperature (or one of the bend- 
ing rigidities e.g. k x ) is varied, we predicted our model 
would follow a generic path like Pi in Fig^. And, 
indeed, these simulationsa observed two specific heat 
bumps, corresponding to two distinct continuous transi- 
tions crumpled-to-tubule and tubule-to-flat [rounded by 
finite membrane size), just as we predictedu. Further- 
more, the shape of the membrane in the phase between 
these two transitions was exactly that of the tubule above 
(see Figure^), and had, within numerical errors, precisely 
the scaling properties and exponents that we predicted 
for phantom tubulesH. Here we present our detailed study 
of these transitions and the tubule phase, in the presence 
of both thermal fluctuations and self-avoidance. 

There are a number of possible experimental realiza- 
tions of anisotropic membranes. One is polymerized 
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FIG. 2. Schematic picture of the tubule phase of 
anisotropic polymerized membrane, with the definition of its 
thickness Rg an d ro ugh ness firms, our predictions for which 
are given in Eqs 



1.1 



and 1.: 



The scaling functions SR t h(x) have the limiting forms: 

x C-Vp/x f or x _> 



Sr(x) oc 



constant, for . 



. constant, 

s h (x) cx <; 



for x 
for x — 



■+ 

00 



(1.4) 



(1.5) 



where v p is the radius of gyration exponent of a coiled lin- 
ear polymer w 3/5. These scaling functions are universal 
(i.e., independent of material parameters and tempera- 
ture), up to an overall non- universal multiplicative fac- 
tor, which can, and will, depend on material parameters 
and temperature. 
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The scaling forms, Eq. 1.4 and [L| imply that for a 
"roughly square" membrane - that is, one with L±_ ~ 
L y = L - in the limit L — > oo 

i? G (L ± ~ L„ = L) oc L" , (1.6) 
^ ms (L ± ~ L y = L) oc L 1 -"^/ 2 , (1.7) 

where we have used the fact that for L y ~ Lj_, the argu- 
ment a; = L y jL\ of the scaling f unct ions Sp^ix) goes to 



infinity as L — > oo, and used Eq.1.3 to simplify Eq.1.7. 

Detailed renormalization group calculations show that 
r] K is strictly positive. Hence, h rms << L for a roughly 
square membrane as L — > oo. Thus, the end-to-end ori- 
entational fluctuations 9 ~ h rms /L oc L~ VkZ / 2 — > as 
L — > oo for such a roughly square membrane, proving 
that tubule order (which requires orientational persis- 
tence in the extended direction) is stable against undu- 
lations of the tubule embedded in d = 3 dimensions. 

On the other hand, in the limit L y >> L±, in which 
the tubule looks more and more like a linear polymer (a 
ribbon of width Lj_ and length L y ), we find 



.3/2 



" rmS ^ r»(3/2-C) 

J-J I 



L 



T 3/2 
l/2+r, K z/2 - U V 



Li, 



L P (L±) 



1/2 



(1.8) 



acting like a rigid polymer with a polymer bending rigid- 
ity 



k p (L±) oc i 



1+zVk 



(1.9) 



It is well knowntfl, of course, that a linear polymer 
does not have long-ranged orientational order i.e., it has 
a finite orientational persistence length Lp. For length 
smaller than Lp(L±) we recover the well-known0 U'J 2 
growth of transverse fluctuations. By equating h rms from 
with the length L y of the tubule, and defining 



1. 



Eq 

(ribbon width-dependent persistent length) Lp(Lj_) to 
be the value of L y at which this equality occurs, we ob- 
tain an estimate for the orientational persistence length 
Lp of a long, skinny tubule: 



Lp(L±) oc L 



1+VitZ 



(1.10) 



We see that only very long, skinny membranes (L y > > 
L±) will be orientationally disordered; for any membrane 
with a reasonable aspect ratio (i.e., L y ~ L±), L y is 
much less than Lp(L±), and the orientational order of 
the tubule persists throughout it. This proves that the 
tubule phase is stable in the thermodynamic limit against 
thermal fluc tua tions. 

Equation [l^ indicates that the effective polymer bend 
modulus k p (L±) is "anomalous", by which we mean the 
fact that k p (L±), grows as a power of L± greater (by 
the "anomalous dimension" T) K z) than 1 (naively ex- 
pected based on dimensional analysis). This together 
with the concomitant an omal ous dimension of the persis- 
tent length Lp(L±), Eq.1.10 embodies the phenomenon 



known as "anomalous elasticity" Jij'0£H9 In addition to 
fluctuating membranes, they have consequences for poly- 
mers whose internal structure is that of a long ribbon of 
dimension Lj_ x L y , with L y » Lp(L±) » L±. Pro- 
vided that L± is large enough that the anomalous elastic- 
ity can manifest itself, the radius of gyration 1Z V G of this 
polymer (which, since L y » Lp, will be coiled) will, 
in fact, grow more rapidly with the transverse dimension 
L± of the polymer than the conventional elastic theory 
would predict. Specifically, we expect: 



Til 



L P {L A 



« L v PL ± 



L P (L ± ) 



1.11) 



while conventional elastic theory would imply TZq oc 

In addition to this anomalous elasticity in the effec- 
tive polymer bend modulus, the fluctuating tubule also 
displays anomalous elasticity for stretching the tubule. 
In particular, experiments that attempt to measure the 
stretching modulus g y of the tubule (def ined mo re p re- 
cisely by the renormalized version of Eqs.^5 and 5.7 ) at 
wavevector q will produce results that depend strongly 
on q, even in the limit q — > 0. In particular, this ap- 
parent wavevector-dependent stretching modulus <7j,(q) 
vanishes as |q| — > 0, according to the scaling law 



<??/(q) = l2 u Sg(q v /qi) , 



(1.12) 



where rj u > is another universal exponent, and S g (x) 
another universal scaling function. 

Similarly, the tubule bend modulus k (also defined 
mo re precisely by the renormalized version of Eqs. |5.4| and 
5.6) becomes strongly wavevector dependent as q — > 0, 
but it diverges in that limit: 



«(q) = % nK S K (q y /qi) , 



(1.13) 



with r] K > yet another universal exponent, and S K (x) 
yet another universa l scalin g function. 
The relations Eqs. 



1.1 



1.2 summarize all of the scaling 



properties in terms of the two universal exponents v and 
z (or equivalently r} K ). Clearly, we would like to predict 
their numerical values. There are three distinct cases to 
be considered, as we decrease the embedding dimension 
d in which the D = 2-dimcnsional membrane fluctuates, 
as illustrated in Fig.|| (the generalization to arbitrary D 
is given in the main text). 



regime III 



regime 



regime I 
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FIG. 3. Illustration (in D — 2) of the three regimes of 
embedding dimension d with qualitatively and quantitatively 
different tubule shape scaling properties. Our estimates of 
d, ~ 6.5 place the physical tubule (d = 3) deep in regime III; 
the strict bound d, > 7/2 guarantees this. 

Regime I: 

For a phantom membrane, or for a membrane with intrin- 
sic dimension D = 2 embedded in a space of dimension 
d > d uc — 11, self-avoidance effects can be asymptoti- 
cally ignored in the tubule phase, and we predicttl 



v 



1 

4 
1 

Z= 2 
Vk = J 
»7« = 1 > 



(1.14) 

(1.15) 

(1.16) 
(1.17) 



Regime II: 

For a self-avoiding membrane with d* < d < d uc = 11 
(with d* > 7/2), we have shown (as we describe in detail 
in Sec. VI) that the bending elasticity is not anomalous, 
i.e., r\ K = 0, as g uaran teed by an exact "tu bule-gauge" 
symmetry (see Sec. VI B ). This, using Eq. 1.3 immediately 
leads to the the exponent relation, z = (1 + 2i/)/3, which 
states that for d > d* all properties of a self-avoiding 
tubule can be expressed in terms of a single radius of 
gyration exponent v. In this range d* < d < d uc = 11 
of embedding dimensionality, the exponents v and z can 
be computed in an e = 11 — d-cxpansion. This has been 
done recently by Bowick and Guitter (BG)H who have 
verified the validity of the Ward identity z = (1 + 2i/)/3 
(for D = 2) perturbatively, to all orders in e. Further- 
more, for all embedding dimensions d > d*, the ab- 
sence of anomalous bend elasticity (i.e., i] K — 0) ren- 
ders the self-avoiding interaction ineffective in stabiliz- 
ing wild transver se t ubule und ulations and for a square 
membranes, Eqs.1.7 and 1.10 show that the D = 2- 



dimensional tubule phase is only marginally stable. For 
D = 2, this d* < d < d uc = 11 regime has: 



2 1 

5 >V> 4> 
z= |(1 + 2i>) , 
'7k = , 

7?u = 3 , 



(1.18) 

(1.19) 
(1.20) 
(1.21) 



Regime III: 



Finally, as we describe in Sec. VI, the physics of the phys- 
ical tubule (i.e., D = 2-dimensional tubule embedded in 
d = 3 dimensions) is much richer than that for the em- 
bedding dimensions d > d», where "tubule-gauge" sym- 
metry imposes strict nonrenormalization of the tubule 
bending rigidity K. For d < d», because of the presence 
of additional elastic nonlinearities (which are irrelevant 
for d near d uc — 11, but become strongly relevant for 



physical dimensionality d < d*), this e-expansion about 
d = d uc = 11 gives no information about the simultane- 
ous role that the self-avoidance and elastic nonlinearity 
play in the physical tubule (D — 2, d = 3 < d*(-D = 2)), 
where they are both important. We find that, as the em- 
bedding dimension d is lowered below d* < d uc = 11 
(d*(D = 2) > 7/2), the nonlinear elasticity becomes 
relevant, destabilizing the fixed point studied in Ref. ||, 
and leading to the breakdown of the z = (1 + 2v)/3 
relation (with the amount of breakdown described by 
a new anomalous elasticity exponent rj K ). Hence phys- 
ical tubules (D — 2, d — 3) are described by a new 
infra-red stable fixed point, that is non-perturbative in 
e = 11 — d, which incorporates the simultaneous effects of 
self- avoidance and nonlinear anomalous elasticity. This 
new fixed point characterizes the d < d* regime (appro- 
priate to a physical tubule) with shape scaling exponents 



v > 



v > 



1 



d- 1 
1 



3 - rj K 

2?7k + Vu = 3 - - 

z 



(1 + 2v) 



(1.22) 
(1.23) 
(1.24) 

(1.25) 
(1.26) 



We cannot calculate exactly the critical embedding 
dimension d*(D) that separates regime II and regime 
III, but we can derive a rigorous lower bound on it 
d*(2) > 7/2. Thus the physical tubule, D = 2, d = 3 
falls in regime III. Our best estimate of d*(2) is that it 
lies between 5 and 7. 

It should be emphasized that all of the exponents are 
universal in a given embedding dimension d. Indeed, for 
d* < d < 11, where all of the exponents are determined 
by the single unknown exponent ^, there are two different 
analytical approximations to v that agree to better than 
1% for d > 8, and to better than 10% for d's greater 
than the likely values of d*. These analytical methods 
are: Flory theorya, which predicts 



v F = 



and the leading order in e = 11 
and GuitterEI, which gives, 



(1.27) 



d expansion of Bowick 



with 



c = 0.13125 



(1.28) 



(1.29) 



We suspect, based on the experience of comparing poly- 
mer exponents obtained from Flory theory with those 
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obtained from the e-expansion, that, although BG's re- 
sults are certainly more accurate near d = 11, when the 
BG and Flory results start to disagree appreciably (i.e., 
below d = 7), the Flory result is probably the more accu- 
rate. Nonetheless, the extremely close agreement between 
these two very different approaches in these high embed- 
ding dimensions increases our faith in both of them. 

In fact, as we describe in detail in Sec. VI B, for D = 2- 
dimensional membrane, is determined by the condi- 
tion that v(d) — > 2/5 as d — ► d+. Using the Flory result 
(Eq. |l.27[ ), this gives d* = 13/2 = 6.5; while using the BG 
result (Eq|L2|) gives d* = 11 - 2/(3c) = 5.92. 



T tubule 
V GV 



7 



3d -5 



(1.31) 



All of the exponents jump discontinuously (as a func- 
tion of d) at d*; figure ^| shows such a plot, schematically, 
for v(d) and rj K {d). 



For a physical tubule, Flory theory, Eq.1.27, implies 



v F (D = 2,d= 3) = 3/4 



(1.30) 



in contrast to the BG result Eq.1.28, which implies 
v t (D = 2,d = 3) = 0.517. What is the correct value 
of v in d = 3? As discussed above, our experience with, 
polymers suggests that Flory theory is more reliabletS 
than the e-expansion when both are pushed well below 
the upper critical dimension. One might be concerned 
that this ceases to be true for tubules, due to the discon- 
tinuous behavior of all of the exponents at , but we will 
present arguments later that suggest that this is not the 
case, and that Flory theory is probably quite accurate in 
the physical Gase-pf d = 3. nn 

It is widelyllZrllS, though not universallyJla^E3 believed 
that self-avoidance destroys the crumpled phase. What 
is definitely known is that the crumpled phase has only 
been seen in simulations of phantom membranes and 
in more recent simulations by Baumgartnenia of a self- 
avoiding plaquette membrane model. It is therefore rea- 
sonable to ask whether our tubule phase will suffer the 
same fate. We think not, for the following reasons: 

1. It is clear that self-avoidance, though a relevant 
perturbation (in physical embedding dimension d < 
d uc = 11) has far less effect on the tubule than 
the crumpled phase, since points on the membrane 
widely separated in the y-direction never bump into 
each other in the tubule phase, but do in the crum- 
pled phase. 

2. The analytic argument that self-avoidance destroys 
the crumpled phase is based pu the Gaussian vari- 
ational (GV) approximatiorOl23, which predicts 

,i . ,i j • r • , crumpled 

that the radius ot gyration exponent v GV — 
4/d, which implies that v > 1 for d < 4, and hence 
that the membrane is extended (i.e. flat) for those 
dimensions (which, of course, include the physical 
case of d = 3). We find that the same Gaussian 
variational approximation leads to the same con- 
clusion for the tubule phase. Our result for D = 2 
is 



and implies v t Qy ule > 1 for d < 4, and hence an in- 
stability of the self-avoiding tubule to an extended 
(i.e. flat) membrane in physical dimensions. 

We are not, however, overly concerned by this re- 
sult, for a number of reasons: 

(a) The Gaussian variational approximation is 
known to be far from trustworthy. For ex- 
ample, it predicts v = 2/d for linear poly- 
mers, which not only is less accurate for all 
d between 1 and 4 than the Flory result 
v = 3/(d + 2), but also incorrectly predicts 
that the lower critical dimension di c below 
which linear polymers are always extended is 
di c = 2, whereas, in fact, it is known exactly 
that di c = 1, a result that is also predicted 
exactly by the Flory theory. Thus, the Gaus- 
sian variational approximation is very unreli- 
able in predicting the lower critical dimension 
of a crumpled object. 

(b) There is a good reason to believe it is equally 
unreliable for our problem as well. If we 
compare the Flory prediction for v with the 
e-expansion calculation of Ref. ^ (which is 
asymptotically exact in d —* 11), in, e.g., 
d = 8, we find they differ less than 1/3 of 1%: 
v e = 0.332 El, VFiory — 1/3 □; while the Gaus- 
sian variational result v l ^ ule = 7/19 = 0.3684 
is nearly 40 times as far off v c as the Flory re- 
sult. This strongly suggests that both Flory 
theory and the e-expansion are more reliable 
than the Gaussian variational approximation, 
and both of them pxedict v substantially < 1 
in d = 3: v F = 3/4 H, v € = 0.517 B. 

(c) Finally, on more general grounds, while the 
Gaussian variational method can be quite use- 
ful, only some of its results can be trusted. 
Certainly it is likely that the trends of, e.g., ex- 
ponents with dimensionality d and D, are cap- 
tured correctly by this theory. The very exis- 
tence of the crumpled phase relies on the pre- 
cise value of v{d) (it disappears if the d < d\ c , 
with di c defined by v{di c ) = 1). However, 
as with any approximate method, especially 
with uncontrolled approximations such as the 
Gaussian method, there is little credibility in 
the actual values of the exponents. Further- 
more, the Gaussian variational approximation 
is very closely related to a large expansion in 
1/d about the embedding dimension d — > oo 
limit £3 It is therefore intrinsically untrustwor- 
thy and ad hoc for small values of d at which 
one is assessing the stability of the tubule (or 
crumpled) phase, which very delicately and 
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sensitively depends on the precise value of v 
at small d. 

In the remainder of this paper we present the details 
of our calculations. In Sec.Q we introduce the Landau- 
Ginzburg- Wilson free energy for our generalized model 
of anisotropic polymerized membranes. In Sec. Ill we 
will first solve this model in mean field theory. From 
this solution we obtain the phase diagram for anisotropic 
polymerized membranes, and identify and characterize 
the new tubule phase as well as the previously studied 
crumpled and flat phases. In Secjl^ we show that the 
scaling properties of the flat and crumpled phases are 
unaffected by the anisotropy. In Secs.|v| and VI we then 
consider the effects of both thermal fluctuations and self- 
avoidance on the tubule phase. We treat this problem 
using Flory theory, renormalization group and Gaussian 
variational methods. We calculate the upper critical em- 
bedding and intrinsic dimensions for both effects, and 
thereby show that both are relevant for the physical case 
of two-dimensional membranes embedded in three dimen- 
sions. We also show that, although there is no anomalous 
elasticity for the bend modulus k along the tubule near 
d = due = 11 (due to aforementioned "tubule gauge" 
symmetry), such anomaly must set in for embedding di- 
mensions d < d*, with d* > 7/2. When this happens, 
the fixed point (perturbative around d = lip which de- 
scribes a self- avoiding (bend elastically raon-anomalous) 
tubule, becomes unstable, and a new fixed point con trols 
the tubul e phase. We derive new exact relations Eqs. |6.62| 
and 6.63 between v and z, which involve anomalous elas- 
ticity exponent r/ K (or r/ u , related to it) and are appro- 
priate for a physical (with anharmonic elasticity) tubule, 
described by this new fixed point. We then use the Florytl 
and extrapolated e = 11 — d-expansiona results for v in 
this relation to determine z and all other tubule shape 
exponents in terms of two constants that, unfortunately, 
we were not able to compute accurately. 

In Se ctio n [y] we also derive the scaling results Eq jl.l| 
and Eq. [l.2| for Rq and h rms , and for the anomalous elas- 
tic theory as well . 

In Section VII we use the renormalization group to 



analyze the crumpled-to-tubule transition. We then con- 
struct a scaling theory of the crumpled-to-tubule and 
tubule-to-flat transitions, and compute within Flory the- 
ory the critic al exp onents for these transitions. 

In Section VIII we summarize, conclude, and make 



some suggestions for further analytic, numerical, and ex- 
perimental work. 



II. MODEL 

Our model for anisotropic membranes is a generaliza- 
tion of the isotropic model considered in Ref. |24. As 
there, we characterize the configuration of the membrane 
by giving the position r(x) in the ci-dimensional embed- 
ding space, of the point in the membrane labeled by a D- 



dimensional internal co-ordinate x. In the physical case, 
d = 3 and D = 2, of course. Throughout the remainder 
of this paper, we will distinguish between D-dimensional 
"intrinsic" vectors and <i-dimensional "extrinsic" vectors 
by using boldface type for the former, and vector arrows 
over the latter. 

We now construct the Landau- Ginzburg- Wilson free 
energy F for this system, by expanding F to leading 
order in powers of r(x) and its gradients with respect 
to internal space x, keeping only those terms consistent 
with the symmetries of the problem. These symmetries 
are global translation invariance r(x) — * r(x) + r D , and 

global rotational invariance r(x) — > M ■ r(x), where f a 

and M are a constant (i.e. x-independent) vector and 
a constant rotation matrix, respectively. Global trans- 
lational invariance requires that F be expanded only in 
powers of gradients with respect to x. We will further- 
more take the membrane to be isotropic in the D — 1 
membrane directions (hereafter denoted by xj_) orthog- 
onal to one special direction (which we call y) . Since the 
physical case is D = 2, this specialization is innocuous. 

The most general model consistent with all of these 
symmetries, neglecting irrelevant terms, is, 



m*)\ = \J d°- l x^dy 



+ K ±y d 2 y r- d 2 ± r + t ± (d^r) 2 + t y (d y rf 

+ ^ (dir- djr) 2 + U -f (d y f. d y rf + u Xy (d^r- d y r)' 

+ ^ Kr- 3 Q V) 2 + v ±y [dirf (d y rf 



+ - I d D x I d D x' <5 (d) (r(x) - r(x')) , 



(2.1) 



where the re's, i's, it's, u's are elastic constants. The first 
three terms in F (the k terms) represent the anisotropic 
bending energy of the membrane. The elastic constants 
t± and t y are the most strongly temperature dependent 
parameters in the model, changing sign from large, posi- 
tive values at high temperatures to negative values at low 
temperatures. Their positivity at high temperatures re- 
flects the membrane's entropic preference for crumpling. 
To see this, note that this crumpled state is one in which 
all the particles in the membrane attempt to cram them- 
selves into the same point r; in this state, the gradients 
with respect to the internal space d^r and d y r seek to 
minimize themselves, which is clearly favorable when t±, 
t y > 0. However, when either of these becomes negative, 
it becomes favorable for the membrane to flatten (i.e., 
extend) in the associated direction, as we shall show in 
a moment. The u and v quartic terms are higher or- 
der elastic constants needed to stabilize the membrane 
when one or both of the first order elastic constants t± , 
t y become negative. Stability requires that 



G 



4 J__L 



> 

> 



> - 



where 



u' xx = v xx + u_lx/(D - 1) 



(2.2) 
(2.3) 

(2.4) 



(2.5) 



The final, b term in Eq.^J represents the self-avoidance 
of the membranes; i.e., its steric or excluded volume in- 
teraction. 



Equation 2.1 reduces to the model for isotropic mem- 

~~ | when t± = t y , 



branes considered in Rcf. 



K±y 



-- o, 

V±y 



h yy ' 
- Av. 



4(6 



K±J_ Ky, 

= 4tt, and 



III. MEAN FIELD THEORY 

We begin our analysis of this model by obtaining its 
mean field phase diagram, at first neglecting the self- 
avoidance interaction. Later, we will consider both the 
effects of fluctuations and self-avoidance. 

In mean-field theory, we seek a c onfiguration r(x) 
that minimizes the free energy Eq. [2.l| (without the self- 

2 

avoidance term). The curvature energies k± (d x r) and 

K y (dyf*) 2 are clearly minimized when r(x) is linear in x. 
We will therefore seek minima of F of the form 



r(x) = (Cxxx,Cy,0,0,...,0) 



(3.1) 



Obviously, uniform rotations r(x) — > M ■ f(x), of any 

such minimum, with M a constant rotation matrix, will 
also be minima. A continuous degenerate set of minima 
is thereby obtained, as usual for a system with a broken 
continuous symmetry. Uniform translations of the entire 
membrane are also allowed, of course. 



Inserting Eq. 3.1 into Eq.2.1, and for now neglecting the 
self- avoidance term, we obtain the mean-field free energy 
for anisotropic membranes 



f = Il 



D-l 



Li, 



\u' xx {D l) 2 Cl + \u yy $ + v Xy (D l)ClC, 2 



, (3.2) 



where L± and L y are the linear dimensions of the flat- 
tened membrane in the _L and y directions, respectively. 

This mean field theory is precisely that studied long 
ago by Fisher et al.c3 for a completely different (mag- 
netic) problem. Minimizing the free energy over and 
( y yields two possible phase diagram topologies, depend- 
ing on whether u'j 



For u'. 



'mi 



xx 
> v 2 



hv > V Xy or U XX* 



ryy < V ±y 



Fig.|l|. Both Cx and ( y vanish for t± 



we obtain the phase diagram in 
This is the 



crumpled phase: the entire membrane, in mean-field the- 
ory, collapses into the origin, (± = = i.e., r(x) = 
for all x. 

In the regime between the positive t±-axis (i.e., the 
locus t y — and t± > 0) and the t y < part of the 
t y = (u yy /v±y)t± line, lies our new y-tubule phase, char- 
acterized by C-L = and Cy — \/\ty\/ u yy > ; t ne mem- 
brane is extended in the y-direction but crumpled in all 
D — 1 ^-directions. 

The _L-tubule phase is the analogous phase with 
the y and _L directions reversed, C, y — and (± = 
■\/\t±\/uj_± > (obviously a symmetrical reversal for 
the physical case of D — 2), and lies between the t± < 
segment of the line t y — (v± y /u' ±± )tj_ and the positive 
t y axis. Finally, the flat phase, characterized by both 

1/2 

C-L = [(\tx\Uyy - \ty\v X y) / '(u' xx U yv - V Xy )] > , (3.3) 
Cy = [{\ty\uxX - \tx\v± y )/{u' xx Uyy - V 2 Xy )] 1/2 > , (3.4) 

lies between the t x < segment of the line t y — 
(u yy /vj_y)tj_ and the ty < segment of the line t y — 

(V±y/U' xx )t±. 

For u' xx u yy < v Xy , the flat phase disappears, and is 
replaced by a direct first-order transition from _L-tubule 
to y-tubule along the locus t y = (v Xy /u' xx )tx (see Fig.^) 



CRUMPLED 




ty > 0. 



FIG. 4. Phase diagram for tethered membranes showing 
our new tubule phase, for the range of elastic parameters 
when the intermediate flat phase disappears. A first-order 
phase transition separates y- and _L-tubule phases. 

The boundaries between the tubule and the crum- 
pled phases remain the positive t y and t± for 
u' xx u yy > v\ y case. 

Note that a direct crumpling transition (i.e. a direct 
transition between the crumpled and flat phase) is very 
non-generic in this picture: only experimental loci that 
pass from t y , t x > through the origin (locus P2 in 
Fig.|l|) can experience such a transition. This transition 
is, in fact, tetra-critical in this picture. 

This does not, however, imply that direct crumpling 
transitions are non-generic. Many membranes will be 
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perfectly isotropic, by virtue of being formed under con- 
ditions of unbroken rotational symmetry (e.g., randomly 
polymerized membranes). As discussed earlier, this set of 
membranes, which is undoubtedly of finite measure, nec- 
essarily lies on the special isotropic subspace of t he full 
parameter space of the model defined by Eq.2^1 speci- 
fied by t± — ty, k±± = Ky, K± y = 0, 
u±± = u± y = Au. and v±± = v± v 



yy = 4 (w + u), 
Av. The val- 



-J-y 



ues of the quartic couplings then satisfy u' ±± u yy > v 
(for u, v > 0), and hence the topology of the phase di- 
agram is Figure [j]. The boundaries of the flat phase for 
those isotropic values of the quartic couplings become 
t y =t±(l + u/v/(D — 1)) and t y = (1 + u/v) t±, respec- 
tively. For u and v both positive (as required by stabil- 
ity), the slopes of these lines are less than and greater 
than 1, respectively; the isotropic locus t y = t±_ there- 
fore lies between the two (i.e., in the flat phase), and 
hence, that model does undergo a direct flat to crumpled 
transition. 

Membranes with any intrinsic bcoken orientational 
symmetry (e.g., in-plane tilt ordeiH, which is quite 
commontS), will generically have t y ^ t±. Furthermore, 
they will not generically have both tj_ and t y vanish at 
the same temperature. A generic locus through the phase 
diagram in Fig.|l| will be like locus Pi, and will necessarily 
have one of the tubule phases intervening between the flat 
and crumpled phases. Our new tubule phase is not only 
generically possible, but actually unavoidable, in mem- 
branes with any type or amount of intrinsic anisotropy. 



IV. FLUCTUATIONS AND SELF-AVOIDANCE IN 
THE FLAT AND CRUMPLED PHASES 



Specifically, k becomes wavevector dependent, and 
«(q) diverges like q~ VK as q — > 0. In .the flat phase 
the standard Lame coefficients /i and A E3 are also in- 
finitely rcnormalized and become wavevector dependent, 
vanishing in the q — ► limit as n(q) ~ X(q) ~ q Vu ; 
the values of rj K and r\ u in the flat phase differ from 
those in the tubule phase, as does their physical in- 
terpretation. The flat phase is furthermore novel in 
that J.L- is characterized by a universal negative Poisson 
raticoEj which for D — 2 is defined as the long wave- 
length limit q -> of a = \(q) / (2(j,(q) + X(q)). The 
transverse undulations in the flat phase, i.e. the mem- 
brane roughness h rms scales with the internal size of the 
membrane as h rms ~ L^, with ( = (4 — D — r] K )/2, 
exactly. Furthermore, an underlying rotational invari- 
ance imposes an exact Ward identity between rj K and ?/ u , 
rj u + 2r/ K = 4 — D, leaving only a single nontrivial in- 
dependent exponent characterizing the properties of the 
flat phase of even anisotropic membranes. The best es- 
timate for rj K in the physical case of a two-dimensional 
membrane (D = 2), embedded in a d = 3-dimensional 
space comes from the self-consistent screening approxi- 
mation (SCSA) of Le Doussal and RadzihovskyEa, who 
find r] K = 4/(1 + vl5) ~ 0.82. The exponent relations 
above then predict r\ u = 0.36 and £ = 0.59. These ex- 
ponents, together with the negative Poisson ratio-predic- 
tions of Le Doussal and Radzihovsky ofcr = — l/3c3 have 
been recently spectacularly verified to high precision in 
very large scale simulations (largest to date) by Falcioni, 
Bowick, Guitter and ThorleifssonEj. 

The root-mean-square (rms) thermal fluctuation 
(|n(x) — z\ ) = (|5n(x)| 2 ) of the local membrane normal 
n(x) about its mean value (here taken to be z) is 



In this section, we show that both the flat and the 
crumpled phases of anisotropic membranes are identi- 
cal in their scaling properties, at sufficiently long length 
scales, to the eponymous phases of isotropic membranes. 

Consider first the flat phase. We can include fluctua- 
tions about the mean-field solutio n by considering small 



deviations from the solution in Eq.3.1 



<|M(x)| 2 ) = <|VMx)| 2 ) , 

= / dW(|Mq)| 2 ) 



oc 



d D q 
n{q)q 2 
oc L 2 -^- D . 



d D q 



(4.2) 



"(x) = (Ci x _l + w_l(x), ( y y + u y {x), h(x)\ , (4.1) where we have imposed an infra-red cutoff q > L 



Inserting this into our initial free energy, Eq.|2Tj with t± 
and t y both in the range in which the flat phase is sta- 
ble, we obtain the uniaxial elastic energy of Ref. ^6|. As 
shown in that reference, fluctuation effects in turn renor- 
malize the anisotropic elastic energy into the isotropic 
membrane elastic energy considered by Refs. |^J^,^. In the 
flat phase, and at sufficiently long scales, the anisotropic 
membranes therefore behave exactly like isotropic mem- 
branes. This in particular implies that the flat phase of 
anisotropic membranes is stable against thermal fluctu- 
ations. As in isotropic membranes, this is due to the 
fact that these very thermal fluctuations dfjMe the bend 
modulus k to infinity at long wavelengthsBQa. 



the integral over wavevectors, L being the smaller of the 
intrinsic linear dimensions Lj_, L y of the flattened mem- 
brane. These fluctuations are finite as L — > oo, when 
2 — r) K — D < 0. In the physical case D = 2, this condition 
is always satisfied since r) R > 0. Thus, membrane orien- 
tational fluctuations remain bounded, and the flat phase 
is stable against thermal fluctuations, for the physical 
case D = 2. Indeed, the SCSA predicts that they remain 
funded down to the lower critical dimension D = \2. 

Note that this stability of the flat phase depends cru- 
cially on the anomalous elasticity, i.e., the divergence of 
K,(q) as q — > 0. In the absence of this effect, which would 
correspond to rj K = 0, the integral over wavevector in 
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Eq. 4.2 would diverge logarithmically for D — 2, describ- 
ing divergent orientational fluctuations leading to an in- 
stability of the flat phase at any non-zero temperature. 
Hence, the flat phase owes its stability to the anomalous 
elasticity (i.e., the fact that rj K > 0). In contrast, as we 
shall show in a moment, the tubule phase is marginally 
stable against thermal fluctuations, even in the absence 
of anomalous elastic effects. Such effects are, nonethe- 
less, actually present for self-avoiding tubules, but they 
are not essential to the stability of the phase. 

Because of this persistent long-ranged orientational or- 
der (i.e., because the membrane is flat), widely intrinsi- 
cally separated parts of the membranes (i.e., points x 
and x', with |x — x'| large) do not bump into each other 
(i.e., never have r(x) = r(x')); hence, the self-avoidance 



interaction in Eq. 2.1 is irrelevant in the flat phase. 

That the crumpled phase of anisotropic membranes is 
identical to that of isotropic membranes is even easier 
to see. When both t± and t y are positive, all of the 
other local terms in Eq. [2.l| , i.e., the k, u, and u-terms, 
are irrelevant at long wavelengths (since they all involve 
more derivatives than the i-terms). Once these irrele- 
vant terms are neglected, a simple change of variables 
x^ = x' \/t_\_/t y makes the remaining energy isotropic. 
Thus, the entire crumpled phase is identical in its scaling 
properties to that of isotropic membranes. 

In particular, the membrane in this phase has a radius 
of gyration Rg(L) which scales with membrane linear 
dimension L like L", with v = (D + 2)/(d + 2) in Flory 
theory, and_vpi^ similar values predicted by e-expansion 
techniques 3 ^ 



V. FLUCTUATIONS IN PHANTOM TUBULES 

In this section, we ignore self-avoidance (i.e., treat 
"phantom" membranes), and consider the effects of fluc- 
tuations on phantom tubules. We will show that these 
fluctuations do not destroy the tubule phase, or change 
the topology of the phase diagram. The detailed prop- 
erties of the tubule phase are, however, modified by the 
fluctuations. 

Let us consider the y-tubule phase (i.e., the tubule 
phase with the tubule axis along the y-axis). To treat 
fluctuations, we perturb around the mean-field solution 
^o(x) = ( y (y, 0) by writing 



^( x ) = {CyV + w(x),/I(x)) , 



(5.1) 



where /i(x) is a d — 1-component vector orthogonal to 
the tubule's axis, which we take to be oriented along the 
y-axis. The average extension factor C, y is near but not 
exactly equal to its mean-field value, because fluctuations 
will change it. Rather, we will choose Cy so that all linear 
terms in ft(x) and it(x) in the resultant elastic free en- 
ergy for these variables are exactly cancelled, in the long 
wavelength limit, by their fluctuation renormalizations. 



This criterion guarantees that /i(x) and u(x) represent 
fluctuations around the true ground state of F. Pre- 
cisely analogous choices have been used in the study of 
bulk smectic A elasticitdLi and the flat-phase elasticity 
of isotropic membrancsOiJ Ll. 



I nser ting the decomposition Eq. 5J into the free energy 
Eq.2.1, neglecting irrelevant terms, and, for the moment 
ignoring the self-avoidance interaction, gives, after some 
algebra, the elastic free energy F tot = F m j t + F e i , where 
F m f t is simply the mean-field free energy for the tubule 
phase 



F m ft - ±L v [tyCy 



(5.2) 



and F e i[u(x),h(x)] is the fluctuating elastic free energy 
part 



F P > = ~ / ,/ 



J-D-l 



x±dy 



+ K{SP y hf + t{d±hf + g^uf 
+ 9y(dyu+^{dyhf + \{d y u) 2 



2-, 



(5.3) 



where k 



t + u± y (y, and 7 = t y 



t = t±+ V ± yCy, 



'■yySy 



gx 



iyy(y are constant coefficients. 
Note first that the coefficient 7 of the linear terms in 
F e i is also the coefficient of the (d y h) 2 term. This is 
a consequence of t he rotation invariance of the origi- 
nal free energy Eq.2.1, which leads to the existence of 
the Goldstone mode d y h. The combination E(u, h) = 
d y u + j{d y h) 2 + \ {d y u) 2 is the only combination of first 

y-derivatives of u and h that is invariant under global 
rotations of the tubule. It is analogous to the nau.- 
linear strain tensor of conventional elasticity theory.EH. 
On these general symmetry grounds, therefore, the free 
energy can only depend on d v u and d y h through pow- 
ers of E{u, h), and this property must be preserved upon 
renormalization. This has two important consequences: 
the first is that, since, as discussed earlier, the coeffi- 
cient of this linear term will be chosen to vanish upon 
renormalization via a judicious choice of the stretching: 
factor £ y , the coefficient of (d y h) 2 will likewise vanish.LJ 
This means that the y-direction becomes a "soft" direc- 
tion for fluctuations of h in the tubule phase. We can 
trace this softness back to the spontaneously broken ro- 
tational symmetry of the tubule state. It is precisely 
analogous to the softness of height fluctuations in the 
flat phase of isotropic membranes, manifested by the ab- 
sence of (d x h) 2 , (dyh) 2 terms in the elastic free energy of 
the flat phase, analogous to Eq. |5.3| (when 7 is tuned to 
0). 

The second important consequence is that the ratios 
of the coefficients of the quadratic (d y u) 2 and the anhar- 
monic d y u(d y h) 2 and (d y h) A terms in F e i must always 
be exactly 4:4: 1, since they must appear together as 







a result of expanding (d y u + \{d y h) 2 + ^(d y u) 2 ) 2 . We 
will show in a few moments that, for this special value of 
these ratios, the long-wavelength anomalous elastic be- 
havior of the "phantom" tubule phase can be calculated 
exactly. 

Recognizing that 7 vanishes after renormalization, we 
can now calculate the propagators (i.e., the harmonic ap- 
proximation to the Fourier transformed correlation func- 
tions) by setting 7 = in Eq.pTa. We thereby obtain 



(/i i (q)/ lj (-q))=fcBT^G ft (q) 
(u(q)u(-q)) = k B TG u (q) , 



where 



G,; 1 (q) = ^iH 

G -l (q) _„.„2 , „ „ 2 



nq y 
9±Q1 + 9v% 



(5.4) 
(5.5) 



(5.6) 
(5.7) 



and bh is a Kroncckcr delta when both indices i and 7 
^ ?/, and is zero if either i or j = y. 

Inspection of the propagators Gh and G u reveals that 
the h- fluctuations are much larger than the u-fluctuations 
for |qj_| ss q 2 , and that it is precisely this regime of 
wavevectors that dominates the fluctuations. Thus, in 
power counting to determine the relevance or irrelevance 
of various operators, we must count each power of |qj_| 
as two powers of q y . It is this power counting that leads 
to the identification of the terms explicitly displayed in 
Eq. |5.3| as the most relevant ones. 

Calculating the root-mean-squared (rms) real space 
positional fluctuations (|/i(x)| 2 ) in the harmonic approx- 
imation by integrating the propagators over all wavevec- 
tors, we find 



OC 



oc 



jD-1 



q±dq y 



1 



q±>L 



(2n) D tq\ 



K <ly 



d D 1 q± T 5/2-D 

q±>L ± 1 q' 



(5.8) 



where we have introduced an infra-red cutoff |q_i_| > LJ 1 
in the last integral. This expression clearly reveals that 
for "phantom" tubules, the upper critical dimension D uc 
for this problem, below which transverse positional fluc- 
tuations diverge is D uc = 5/2; this in pri nciple (but see 
discussion of dominant zero modes in Sec. VI B ) allows a 
quantitatively trustworthy e = D uc — 2~ 1/2 expansion 
for the physical membrane of D = 2. This should be con- 
trasted with the result D uc ™4 for the analogous critical 
dimension in the flat phaselB 

The lower critical dimension D\ c below which the 
tubule is necessarily crumpled in this problem is also low- 
ered by the anisotropy. Considering the fluctuations of 
the membrane normals V/i in the harmonic approxima- 
tion, one sees immediately that the largest of these is the 
fluctuation in the y-direction, 



(|<5n y (x)| 2 > = (\d y h(x)\ 2 ) 



d D x q±dq y q y 



q±>L - ± > (2tt)D tql + Kq y 

\ V ~ " . (5.9) 



a q±. T 3/2-D 



q±>L ± 1 ql 

which clearly only diverges in the infra-red L±_ — ► 00 limit 
for D < Die — 3/2 (but again see discussion of dominant 
zero modes in Sec.lVI Bj). 



In the argot of the membrane field, the elasticity of 
phantom tubules is anomalous. In contrast to the flat 
phase, however, for phantom tubules, the exponents char- 
acterizing the anomalous elasticity can be calculated ex- 
actly. To see this, we first note that the w-fluctuations go 
like 1/q 2 in all directions and hence are negligible (in the 
relevant wavevector regime |qj_| « q 2 ) relative to the h- 
fluctuations which scale like 1/q in this regime. This jus- 
tifies neglecting the ^(d y u) 2 piece of the invariant E(u, h) 
operator. This also emerges from a full renormalization 
group treatments, which shows that this term is strongly 
irrelevant. Once it is neglected, the elastic free energy is 
quadratic in u, and these phonon modes can therefore be 
integrated exactly out of the partition function 



Z 



VuDhe~P F ^ uM . 



(5.10) 



Once this is done, the only remaining anharmonic term in 
the effective elastic free energy for h is, in Fourier space, 

F anh [h] = ~ [ (Ufa) ■ h(k 2 )) (h(k 3 ) ■ h(k 4 ) 

4 Jk u k 3 M ' V 



k yl k y 2k y3 k y 4V h (q) 



(5.11) 



where q = k± + k 2 and ^1+^2 + ^3- 
vertex Va(q) above reduces to 



V h (q) = 



9y9±qi 



g y qy 



9xq A 



-kA = 0. The effective 



(5.12) 



whi ch is irrelevant near the Gaussian fixed point (but see 
Sec. VI B), as can be seen by the simple anisotropic power 
counting described above. 

The exact cancelation of the relevant terms in F an h [h] 
above is a direct consequence of the 4:4:1 ratios of the 
coefficients of the quadratic (d y u) 2 and the anharmonic 

d y u(d y h) 2 and (d y h) A terms in F e i that was discussed 
earlier. Given this cancelation, F an h[h] is now clearly 
less relevant than the anharmonic vertices d y u(d y h) 2 
and (dyh)^ in the original free energy (before we inte- 
grated out the phonons u). This is because the factor 
Vft,(q) oc q\/{g y q 2 + 9±q±) vanishes like q 2 in the rele- 



vant limit |qj_| ~ q 2 , q y — > (the other factors in Eq.5.11 

are precisely the Fourier transform of (d y h)^, of course). 
This lowers the upper critical dimension for anomalous 
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elasticity of the h field to D uc = 3/2. Thus, in the phys- 
ical case D = 2, there is no anomalous elasticity in h ; 
that is, the elastic constants t and k in Eq.5.6 are finite 
and non-zero as q y 



0. 



However, as asserted earlier, the full elasticity Eq. 5.3 



before u is integrated out, is anomalous, because g y is 
driven to zero as q y — > 0. Indeed, a self-consistent one- 
loop perturbative calculation of gyja) , obtained by eval- 
uating the Feynman graph in FigJq, gives 



5y(q) =9y 



(5.13) 



k B Tg 2 y (q) P 2 y (p y - q y ) 2 d D - 1 P± d Py /(2TT) 



D 



to 



4p)p» 



(*Ipj 



•«(|p-q|)(p» 



where g° is the "bare" or unrenormalized value of g y . 




FIG. 5. Feynman graph equation for the self-consistent 
evaluation of <? H (q). 

Our earlier argument shows that k(p) can be replaced 

, since the h elast icity 



by a constant in Eq. 5.13 
is not anomalous 



as p 



The self-consistent equation Eq. 5.13 
can be solved by the ansatz, 



g y (<l) = q^S g (q y /q z ± ) . (5.14) 
Simple power counting!! then shows that we must choose 



Z= 2 



Vu 



B 

i 

5 -2D 



(5.15) 
(5.16) 



It is straightforward to verify that these results hold to 
all orders in perturbation theory; that is, at every order, 
the leading dependence on q of the contribution to g y 

1 /2 

scales like q^"S g (q y /q^ ) with r) u = 5 — 2D. 

It is straightforward to verify to all orders in perturba- 
tion theory that there is no such renormalization of g±. 
This is because of the anisotropic scaling q± ~ q 2 , which 
implies that all vertices proportional to powers of perpen- 
dicular gradients of h, i.e., powers of Vj_/i are irrelevant. 
Since only such vertices can rcnormalizc g±\ V ±u\ 2 , there 
are no relevant renormalization of g±. As a result, g± re- 
mains finite and non-zero, or, in a word, non- anomalous, 
as |q| — > 0. 

Using the facts that g y (q) is independent of q^ as 
|qjj = q± —> for fixed q y , and, likewise, to be inde- 
pendent of q y as q y — > for fixed q± , we can obtain the 
limits of the scaling function S g (x) : 



S g (x) cx 



constant, x 



x 



X 



■ OO 

. 



(5.17) 



For phantom membranes with D 
z = 1/2, so we find: 



Vu 



1 and 



(5.18) 



We will now use this result to compute the mean- 
squared real space fluctuations ((u(L±,y) — u(0,y)) ) = 
(Au 2 ) of u(x). These can be obtained via the equiparti- 
tion theorem and by summing all of the Fourier modes, 
yielding: 



(Au 2 



1 



q±>L ± , q y >L y 



dq±dq y 

L- 1 ( 27r ) 2 9y(<l)qy +9±<li 



(5.19) 



Let us assume, and verify a posteriori, that the inte- 
gral in this expression is do minate d by wavevectors with 
q y « y/g±- Then, using Eq. 5.18, we see that 



(Au 2 



dq±dq y 



1 



?_ L >L- 1 ,g H >L- 1 (2tt) D cjqZql + g ± ql ' 



(5.20) 



where c is a constant. Inspection of this integral reveals 
that it is dominated by q's for which the two terms in the 
denominator balance; this means q y ~ q 3 / 4 << -JqZ, the 
last extreme inequality holding as |q| — > 0. This verifies 
our earlier a posteriori assumption that q y << ^fqZ m 
the dominant wavevector regime. 

Now, changing variables in the integral q± = Q±/L±, 
} y /L 3 / 4 , we find 



1v 



where 



(Au') = L^S^Ly/LT) , 



3/4, 



(5.21) 



S u (x) 



dQ±dQ y 



Q±>h 



(2t:) 2 c^QZQl+ gi 



(5.22) 



We note that the scaling form for the u phonon correla- 
tions is different tha n th at o f th e height field h as sum- 
marized in e.g. Eqs.1.1 and 1.2, and discussed in more 



detail below. 

The limits of S u (x) scaling function can be obtained 
just as we did for S g (x); we find, including "zero modes" 
(see below): 



, constant, x — > oo 
S *( x )<x{ x -i X ^Q. 



(5.23) 



For roughly square membranes, L y ~ L±_ = L, so, as 

L — > oo, L y /L 3 / 4 —> oo, and the first limit of Eq.5.23 is 
the appropriate one. This gives 
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(Am 2 ) oc L 



1/4 



(5.24) 



The authors of Ref. |2j measured a quantity that 
should scale like (Au 2 ) in their simulations of a square 
anisotropic membrane. They did this via their vividly 
named "salami" method: measuring the moment of in- 
ertia tensor of a "salami" slice, a set of N points that 
all had the same internal y coordinate (for a y-tubule 
phase). It is straightforward to show that the smallest 
eigenvalue of this tensor should scale like N(u(x.) 2 ), since, 
as we shall see in a moment, the mean squared displace- 
ments in the other directions are much b igger than those 
in the y direction. Therefore, from Eq.5.24 we predict 
that the smallest eigenvalue of this salami slice moment 
of inertia tensor scales like NL 1 / 4 . BFT actually fit this 
eigenvalue to AT log L, which might appear to disagree 
with our prediction, until one recognizes that for L's be- 
tween 32 and 100 (where most of the data of Ref. || is 
taken), L 1 / 4 = (e/4) \ogL to an accuracy of better than 
1%. Thus, their fit is certainly consistent with our pre- 
dictions. T o mo re s trenu ously test our full scaling pre- 
dictions Eq.5.21 and 5.23, one could simulate membranes 
with aspect ratios quite diffe rent from 1. In particular, 
we predict based on Eq.5.21 that increasing L y at fixed 
Lj_ from an initially square configuration would not in- 
crease this smallest eigenvalue; nor would decreasing L y 

decrease it, until an aspect ratio L y ~ L 3 / 4 is reached, 
beyond which this eigenvalue would increase like L y . 

We now turn to the computation, for the phantom 
tubule, of the tubule radius of gyration Rq and roughness 
hrma, defined by: 

Rl= (\h(L ± ,y)-h(0 ± , y) | 2 ) , (5.25) 
h 2 rms = (\h(x ± ,L y )-h(x x ,0)\ 2 ) , (5.26) 

where L_l spans the intrinsic _L space of the membrane. 
Because Rq is by definition the root-mean-square (rms) 
distance between two points at the same y, it is roughly 
the radius of a typical cross-section of the tubule per- 
pendicular to the tubule axis. Likewise, h rms measures 
fluctuations between points widely separated along the 
tubule axis; hence, it gives the polymer-like transverse 
"wandering" of the tubule. See Fig.|| for an illustration 
of Rq and h rms . 

The reason we distinguish between these two quantities 
is that they scale in different ways with the membrane 
dimensions Lj_ and L y , in contrast to one's naive expec- 
tations. This happens because there are large contribu- 
tions to both quantities from " zero modes" , by which we 
mean Fourier modes with either or q y = 0. Those 
with = correspond to polymer-like undulations of 
the entire tubule. Recognizing the existence of both types 
of modes, we Fourier decompose h(x) as follows: 



ft(x) 



1 



1 Ly 



1 



h 0y (q y )e^y 



-D-l 



fto_L(q_L)e 1 ' 



(5.27) 



where B, Oy, and _L denote "bulk modes" (i.e., modes 
with neither qj_ nor q y = 0), and "zero modes" (i.e., 
modes with either or q y — 0), respectively. Note 
that we have chosen different normalizations for the three 
types of modes. For phantom membranes, we proceed by 
inserting this Fourier decomposition into the harmonic, 
h dependent piece of the elastic free energy F e i (which is 
justified, since, as shown above, the elasticity for h for a 
phantom tubule is not anomalous), obtaining: 

q 



Qy 



(5.28) 



± 1 an d L y 



Note the explicit presence of the factors of L 
for the 0- modes. Applying equipartition to Eq. 5.28 , we 
can obtain the mean squared fluctuations of the Fourier 
modes: 



IMq)| 2 ) 



k B T{d-D) 



tq\ 



nql 



(\ h Qy(q v )\ 



2 v _ k B T{d-D) 



(\h ±(<l±)\ 2 ) = 



k B T(d-D) 
Lytq 2 ± 



(5.29) 



(5.30) 



(5.31) 



Using these expressions inside Eqs.5.25 and 5.2£, and be- 
ing careful about converting sums on q into integrals, we 
get 



R 2 G = 2(d-D) 



k R T 



iD-l, 



\D-X 



-y JL- 1 (2tt 

d D ~ 1 q±dq y (1 — e 



tq± 



(27T) 



D 



tq 



nql 







(5.32) 



where the subscripts L ^ and L y 1 denote infra-red cut- 
offs |qx| > L^ l ,q y > Ly 1 , with L ± = \L ±\. 

We observe here that Rq in Eq. 5.32 does not receive 
any contribution from the q± = "zero mode" (i.e., in 
addition to the bulk mode, Rq receives a contribution 
only from the q y = "zero mode" ) . 

Scaling L± out of both integrals for Rq by the change 
of variables = q±L± and Q y = q y \J L± , we obtain 



,CiL 



3-D 



L h 



lt- d ir(4^)) , 



(5.33) 
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where 



Ci = 2{d-D)k B T 



jD-1 



q± (1 



o»qj_-Lj_ x 



/l (27T 

is a constant of O(l), and 
I R (x) = 2{d-D)k B T 



(5.34) 



d D - 1 q_ L dq y (1 - e'l 111 ) 



(2^ 



(5.35) 



with the unit vector along Lj_. Defining the scaling 
function 



S R {x) = \l^- + I R {x) 



(5.36) 



we see that Rg can be rewritten in the scaling form 



Rg(L±, L y ) — L±S r(L v I Xj_) 
with, for phantom membranes, 
5 -2D 



(5.37) 

(5.38) 
(5.39) 
continues 

to apply when self-avoidance is included, but with dif- 
ferent values of v and z, and a different scaling function 
Sr(x). For phantom membranes, from our explicit ex- 
pression for the scaling function Sr, we see that it has 
the limiting forms: 



We will see later that the scaling form Eq. 5.37 



S R (x) OC 



X 



constant, 



for x 
for x — 



■* 

oo 



(5.40) 



In particular, the limiting form as x — > oo implies that 
for the physically relevant case of a square membrane 



L 



L , 



dominate, and we obtain 



oo, for which L„, 



» L z ± , bulk modes 



R G oc L\ 



(5.41) 



The simulations of BFTa measured Rg for phantom 
tubules by calculating the largest moment of inertia for 
for a set of membrane points that all had the same value 
of the intrinsic coordinate y. While we have used here 



a slightly different definition, Eq.5.25, the square root of 
this moment of inertia should scale like our Rg- And, 
indee d, BF T found that it did scale like a power of L, as 
in Eq. 5.41 , with v = 0.24 ±0.02 in excellent qu antit ative 
agreement with our predictions of v = 1/4, Eq.5.38 eval- 
uated in D = 2. It would be of great int erest to test our 
full anisotropic scaling prediction of Eq. |5.37| by varying 
the aspect ratio of the membrane in such simulations. 
For instance, one could fix L± and increase L y \ we pre- 
dict that one should observe no change in Rq ■ The same 



should hold if one decreased L y at fixed L 



Rg should 

remain unchanged until L y ~ y/L±, at which point the 
tubule should begin to get thinner (i.e. Rg should de- 
crease) . 



Equations 5.36 and 5.40 also correctly recover the limit 
of L y = constant << L\ — > oo, where the q y = "zero 
modes" dominate, the tubule simply becomes a phantom, 
coiled up, D — 1-dimensional polymeric network of size 
L± embedded in d— 1 dimensions, with the radius of gy- 



ration Rg{L±) ~ L 



(3-U)/2 



In the physical dimensions 



(D = 2 and d — 3) in particular this gives a coiled up 



ideal polymer of length L± with Rg - 



L 1 / 2 , as expected. 



We now turn our attention to the calculation of the 
tubule roughness h rms . As we will see, here the = 
zero mode will play an essential role and will dominate 
the transverse undulations for "very long" tubules, which 
(because of anisotropic scaling) in particular includes 
tubules made from square membranes. Using the defi- 



nition of h rms , Eq. 5.26| , we have 



h 2 rms = 2(d - D) 



L 



D-l 



dq y 



Ly 1 ( 27r ) K (ly)<ly 



+ k B T 



jD-1 



qj_dq y (1 



4 
q v L 



(1 



Jq y L y 



tqj 



Kq: 



(5.42) 



Here we observe that h rms in Eq. 5.42 does not receive 
any contribution from the q y — "zero mode" (i.e., in 
addition to the bulk mode, h rms receives a contribution 
only from the qj_ = "zero mode"). This is to be con- 
trasted with the behavior of Rg that we noted following 
Eq.5.32, and is responsible for the differences in scaling 



properties of Rg and h rms , notes above. 

Now, for perverse and twisted reasons of our own, 
we choose to scale L yi rather than L±, out of the in- 
tegrals in this expression, via the change of variables 
Q y = q y L y , Q± = q±L y , which leads to 



h 2 = 



= ( c * L y 



L 



D-l 



5-2Z5 

J y 



h{^=)) 



(5.43) 



where 



C 2 = 2(d-D)k B T f^J— 

Jl 27T K 

is yet another constant of O(l), and 
I h {x) = 2{d-D)k B T 



dQ y (l-ep (5 44) 



d D - 1 Q ± dQ y (l-e* 



2;1 {2n) D tQl + KQ y ' 
(5.45) 

Defining the scaling function 



S h (x) = yJc 2 x* D -V+I h (x), (5.46) 
we see that h rms can be rewritten in the scaling form 
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hrms(L±,L y ) — L^Sh{L y /L z j_) 
with, for phantom membranes, 
, 5 -2D 

1 



(5.47) 

(5.48) 
(5.49) 



Again, this scaling law Eq |5.47 continues to apply when 
self- avoidance is included, but with different values of £ 
and z. 



Equations 5.37 and 5.47 give information about the 
tubule roughness for arbitrarily large size Lj_ and L y , 
and arbitrary aspect ratio. For the physically relevant 
case of a square membrane L\ 
which L y >> we obtain, 



L 



L 



C+(D-l)/2z 



oc 



L 



(D-l)/2 



OC L C + (D-l)(l-z)/2z 



oo, for 

(5.50) 
(5.51) 



Equations 5.48| , 5.49 then give, for a D = 2 phantom 
tubule, C = 1/2, z = 1/2 



7-3/2 

^y 

Irms ~ 1/2 
i 



and therefore predicts for a square membrane 



(5.52) 



(5.53) 

This prediction for square phantom membranes has 
also been spectacularly quantitatively confirmed in sim- 
ulations by BFT0. Their ingenious procedure for deter- 
mining h rms is rather involved, and the interested reader 
is referred to their paper for a clear and complete dis- 
cussion of it. The bottom line, however, is that they 
find h rms ~ L 7 , (our 7 is £ in their notation) with 
7 = 0.895 ± 0.06, in excell ent agreement with our pre- 
diction 7 = 1 from Eq. [5.53 above. As with Rg, it would 
be interesting to test the full scaling law Eq.5.47 by sim- 



ulating non-square membranes, and testing for the in- 
dependent scaling of h rms with L v and L±. Note that, 



unlike Rg, according to Eq.5.52 



diate growth (reduction) when one increases (decreases) 



L y at fixed L 



Because, unlike the flat phase, no log(L/a) correction 
arises, the (D = 2) phantom tubule is just marginally 
stable, but with wild transverse undulations which scale 



linearly with its length. As we will see in Sec. VI, these 
wild fluctuations will be suppressed when the effects of 
self- avoidance are included. 

The above discussion also reveals that our earlier con- 
clusions about the lower critical dimension D\ c for the 
existence of the tubule are strongly dependent on how 
L± and L y go to infinity relative to each other; i.e., on 
the membrane aspect ratio. The earlier conclusion that 



Die = 3/2 only strictly applies when the bulk modes 
dominate the physics, which is the case for a very squat 



membrane, with L y f» L z 



in which case L y << Lj_. 



For the physically more relevant case of a square phan- 
tom membrane, from the discussion above, we find that 
Die = 2~ , where the ~~ superscript means that there are 
no logarithmic corrections at D = 2 and therefore strictly 
speaking the D = 2 tubule is marginally stable. 

Equations 5.37 and 5.47 also correctly recover the limit 



of Lj_ — constant « L y — > 00, where the tubule si mply 



becomes a polymer of thickness Rq(L±) given in Eq.5.25 
of length L y embedded in d — 1 dimensions. As already 
discussed in the Introduction for a more general case of 
a self-avoiding tubule, these equations then correctly re- 
cover this polymer limit giving 



hrms ~ Lp{Ly/ Lpf'l 2 , 

with Lj_-dependent persistent length 
Lp(Lx) oc L^- 1 . 



(5.54) 



(5.55) 



which agrees with Eq.1.10 of SccM for D = 2 when one 



remembers that, for the phantom membranes, r\ K = 0. 
So, as expected for a phantom tubule, if L± does not grow 
fast enough (e.g. remains constant), while L y — > 00, the 
tubule behaves as a linear polymer and crumples along 
its axis and the distinction between the crumpled and 
tubule phases disappears. 

To summarize: the radius of gyration Rq and the 
tubule roughness h rms scale differently with membrane 
size L for a square membrane because the former is dom- 
inated by bulk modes, while the latter is dominated by 
q± = "zero modes" . 



VI. SELF- AVOIDANCE IN THE TUBULE PHASE 

We now look at the effects of self-avoidance on the 
tubule phase, and begin by calculating the upper critical 
embedding dimension d uc below which the self-avoidance 
becomes relevant in the tubule phase. A model of a self- 
avoiding membrane in the the tubule phase is described 
by a free energy functional whi ch is a combination of the 
elastic free energy F e i from Eq.5.3 and the self-avoiding 
interaction Fsa from Eq.2.1 specialized to the tubule ex- 



tended in y-direction using Eq.5.1 for r(x) 



F sa = \J dydy' Vd B - V_l 5^ (h(x ± ,y) - h(x? ± ,y 

x 5(C y y + u(x±,v) - Cyy' - u ( x 'j_,y')) ■ (6- 

If the in-plane fluctuations u scale sub-linearly with y 
(which we will self-consistently verify a posteriori that 
they do), at long length scales one can ignore the phonons 
inside the self-avoiding interaction above. This can be 
confirmed more formally by an explicit renormalization 
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group analysis^. We then obtain a self-avoiding inter- 
action that is local in y, with corrections that are irrel- 
evant in the renormalization group sense and therefore 
subdominant at long length scales. The appropriate free 
energy that describes a self-avoiding tubule is then given 
by 



F = - 



_dy 



d y u- 



\{d v hf 



+ K (& h f + t(d±hf+g ± (d±uf 



9y 



9 v u+-{dyhY 



+ v ldyd D - 1 x ± d D - 1 x' ± d( d - 1 \h(x±,y)-h(x'^y)) , 

(6.2) 

where v = b/2(, y . 

It is important for simulators to note that, although 
the self-avoiding interaction is effectively local in intrinsic 
coordinate y, this does not mean that the effects of self- 
avoidance can be included in simulations that have each 
particle on the membrane avoid only those labeled by the 
same intrinsic y coordinate. Such a simulation, rather, 
models the very different (unphysical) self-avoiding inter- 
action 



dyd l x^d 1 x'j_5(u(x±, y) — u(x'j_, y) 



x ^- 1 '(ft(xj.,y)-ft(x , ±J y)) 



t SA 



(6.3) 

which accounts for interaction only of particles that have 
the same intrinsic coordinate y and the same extrinsic 
coordinate. For large membranes, this unphysical inter- 
action is smaller than the true self-avoiding interaction in 
Eq. |6.2| by a factor that scales like the inverse of the rms 
fluctuations of u: (u 2 )" 1 / 2 , as can be se en t rivially from 
the scaling of the 5-function of u in Eq.6.3. Since these 
fluctuations of u diverge as L± — > oo like u rms ~ L ™ , 
with Cu > (e.g., C« = 1/8, for d > 11 and D = 2), the 
wrong self-avoiding interaction in Eq. |6.3| drastically un- 
derestimates the true self-avoiding interaction by a factor 
that diverges in the thermodynamic limit. Although it 
is tempting to do so in simulations, one must be careful 
not to implement the unphysical self-avoiding interaction 
in Eq |6.3| . Since it might be difficult to implement the 
approximate (but asymptotically exact) self-avoiding in- 
teraction of Eq. |6~2| in simulations, it is easiest to simulate 
the unapproximated interaction in Eq. |6.l| . 

In the next three subsections we analyze the properties 
of a self-avoiding tubule described Jay this nonlinear elas- 
tic free energy, using Flory theorya, the renormalization 
group, and the Gaussian variational methodtl. 



A. Flory theory 

The effects of self-avoidance in the tubule phase can 
be estimated by generalizing standard Flory arguments 



from polymer physics^ to the extended tubule geometry. 
The total self-avoidance energy scales as 



where 



E SA cx Vp 2 



V cx R d Q X L v 



(6.4) 



(6.5) 



is the volume in the embedding space occupied by the 
tubule and p = M/V is the embedding space density of 
the tubule. Using the fact that the tubule mass M scales 



like L 



D- 



1 L y , we see that 



E S a oc 



L n ,L 



2(D-1) 



r>d-l 
G 



(6.6) 



Using the radius of gyration Rq cx L v ^ and consider- 
ing, as required by the anisotropic scaling, a membrane 
with L± oc L y , we find that E$a oc Ly SA around the 
phantom fixed point, with 



\ SA = l+4(D-l)-2(d-l)v 



(6.7) 



Self-avoidance is relevant when \$a > 0, which, from 
the above equa tion, happens for v — v p h = (5 — 2D) /A 
(as per Eq.p.3q) when the embedding dimension 



d < df* 



QD - 1 
5 -2D 



(6.8) 



For D = 2-dimensional membranes, df^ = 11. Thus, 
self- avoidance is strongly relevant for the tubule phase 
in d = 3, in contrast to the flat phase. 

We can estimate the effect of the self-avoidance inter- 
actions on Rq{L i ) in Flory theory, by balancing the es- 
timate Eq. p.6| for the self-avoidance energy with a similar 
estimate for the elastic energy: 



E, 



last/ 



Rg 

Li 



L i L, 



(6.9) 



Equating E e i ast i C with Esa-i we obtain a Flory estimate 
for the radius of gyration Rq: 



R G (L±_) cx L v l , v F = 



D + l 
d + 1 



(6.10) 



which should be contrasted with the Flory estimate of 
Vp = (D + 2)/(d + 2) for the crumpled phase. The sim- 
ilarity of the expressions is not surprising, since for the 
tubule phase the y-dimension decouples in both the in- 
trinsic and the embedding spaces and is not affected by 
the self- avoidance. For the physical case D = 2, d = 3 



Eq.6.10 gives 



Rg cx Lj_' 



(6.11) 
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a result that is known to be exact for the radius of 
gyration ofp-a D = 1-polymer embedded in d = 2- 
dimensions.E3 Since the cross-section of the D = 2- 
tubule, crudely speaking, traces out a crumpled polymer 
embedded in two dimensions (see Fig.^J), it is intriguing 
to conjecture that v = 3/4 is also the exact result for 
the scaling of the thickness of the tubule. Unfortunately 
we have no strong arguments supporting this appealing 
conjecture. 

For a square membrane, L y ~ L±, it is straightforward 
to argue, as we did previously, that the q y = zero modes 
do not co ntrib ute to Rg, and Lj_ is the relevant cutoff. 
Hence Eq.6.10 gives the correct radius of gyration. More 
generally, we expect 



R ( ,(L__. L„) x r>; ( ^ 



(6.12) 



where Sr(x) is the scaling function give n in Eq. |l.5| and 
z is the anisotropy exponent given in Eq.1.3. 



B. Renormalization group and scaling relations 

In this subsection, we present a renormalization group 
analysis of the physical self-avoiding membrane, which 
will also require a simultaneous treatment of the non- 
linear elasticity that was already present in a phantom 
membrane, as discussed in Sec.|v|. 

The correct model, which incorporates the effects of 
both the self-avoiding interaction and th e a nharmonic 
elasticity, is defined by the free energy Eq. 6^ . 



F = - I d^xydy 



n{dlhf+t{dihf+g^uf 



+ g y [d y u+ -{d v h)' 



+ v jdyd D - 1 x ± d D - 1 x' x S^(h{x ± ,y)-h(x' x> y)) , 

(6.13) 

where we have set 7 = and dropped the subdominant 
phonon anharmonicity. 

It is convenient for the purposes of this section to 
choose the units of length such that t = k = 1 through- 
out, and choose the renormalization group rescalings to 
keep them fixed at 1 even after the diagrammatic cor- 
rections are taken into account (i.e., beyond the tree- 
level). We, follow the standard renormalization group 
procedureLJ: 

(i) Integrate out fluctuations of the Fourier modes w(q) 
and h(q) of the fields it(x) and h(x.) with wavevectors in 
the high wavevector shell Ae~ l < q± < A, — 00 < q y < 
cxd, where the ultraviolet cutoff A is of order an inverse 
microscopic length, and I is a parameter known as the 
"renormalization group time" . This integration can, of 



course only be accomplished perturbatively in the non- 
linear couplings v and g y . 

(ii) Anisotropically rescale lengths (xj_, y) and fields 
(/i(x), u(x)), so as to restore the ultraviolet cutoff to A: 

(6.14a) 
(6.14b) 
(6.14c) 
(6.14d) 



= e z x' 



e zl y' 



h(x) = e vl h'(x') , 
u(x) = e (2,y ~ z V(x') 



where we have chosen the convenient (but not necessary) 
rescaling of the phonon field u so as to preserve the form 
of the rotation-invariant operator (d y u + ^(d y h) 2 ) 2 . 

(iii) Define the effective length-scale dependent cou- 
pling constants so as to bring the resulting long wave- 
length effective free energy into the same form as Eq. |6.13 



As discussed above, we will choose the arbitrary rescal- 
ing exponents v and z so as to keep the renormalized k(1) 
and t(l) equal to one. This choice of v and z can be shown 
by standard renormalization group arguments t o be the 
v a nd z that appear in the scaling function Eqs.1.1 and 
1.2, as we will demonstrate later in this subsection. 

The result of the three steps of the above renormaliza- 
tion group transformation (i.e., mode integration, rescal- 
ing, and coupling redefinition) can be summarized in dif- 
ferential recursion relations for the flowing coupling con- 
stants: 

dt 

dl 
dn 
~dJ 
dg y 
dl 

dg± 
dl 
dv 

~dl 

where the various /-functions represent the graphical 
(i.e., perturbative) corrections. Since the self-avoiding 
interaction only involves h, and the parameters in the 
h propagator (t and k) are going to be held fixed at 1, 
the graphical corrections coming from self-avoiding inter- 
action alone depend only on the strength v of the self- 
avoiding interaction. Therefore, to all orders in v, and 
leading order in g y , ft{v) and f v (v) are only functions of 
v and f K (g y ,g±) and f g {g y ) are only functions of g y and 
9±- 

It is important to note that g± suffers no graphical 
corrections, i.e., Eq.6.18 is exact. This is enforced by an 
exact symmetry 



[2v + 


z- 


-D- 


3- 


ft(v)]t , 


(6.15) 


[2v- 


3z 


+ D 


- 1 - 


r- fn{9y,9x)]K , 


(6.16) 


[4i/- 


3z 


+ D 


- 1 - 


- fg(9y)]9y , 


(6.17) 


[4i/- 


z- 


D- 


3k 


L ) 


(6.18) 


[2D - 


- 2 


+ z- 


(d- 


- l)f - fv(v)]v , 


(6.19) 



«(x_l, y) -> u(x_l, y) + x( x ±) , 



(6.20) 



where x( x J-) is & n arbitrary function of xj_, under which 
the nonlinearities in F are invariant. 

We further note that there is an additional tubule 
"gauge" -like symmetry for g y = 



16 



h{-x._L,y) -> h(pc±,y) + <j>(y) 



(6.21) 



under which the only remaining nonlinearity, the self- 
avoiding interaction, being local in y, is invariant. 
This "tubule gauge" symmetry demands that f K {gy = 
0,g±) = 0, which implies that if g y = 0, there is no 
divergent renormalization of n, exactly, i.e., the self- 
avoiding interaction alone cannot renormalize k. This 
non-renormalization of k by the self-avoiding interac- 
tion, in a truncated (unphysical) membrane model with 
g y = 0, has been recently verified to all orders in a per- 
turbative renormalization group calculation^. 

To see th at the v a nd z obtained as fixed point solu- 
tions of Eqs. 6.15 - |6.19 have the same physical significa nce 
as the v and z defined in the scaling expressions Eqs jl.l| 
and |l.2j for the radius of gyration Rq and tubule wig- 
glyness h rms , we use the renormalization group transfor- 
mation to relate these quantities in the unrenormalized 
system to those in the renormalized one. This gives, for 
instance, for the radius of gyration 

R G (L ± ,L y ;t(0),K(0),...) = 

= {\K(L±,y)-h(p±,y)\ 2 ) 1 ' 2 

L„,t(0),re(0),... 

= e» l (\h(e- l L ± ,y)-h(0 ± ,y)\ 2 ) 1 / 2 



e- zl L y ,t(l), K (l),... 



e vl R G (e- l L x ,e-* l L y ;t(l), K (l),...) 



(6.22) 



where t(l), k(1), . . . stand for all flowing coupling con- 
stants whose evolu tion with / is determined by the recur- 
sion relations Eqs, 6.15| - 3.19 . Choosing I = U = logLj_ 
this becomes: 



Rg(L±,L v ; t, k, . 



LlR G (l,L y /Ll;t(Q,K(Q,. 



(6.23) 



This relation holds for any choice of the (after all, arbi- 
trary) rescaling exponents v and z. However, if we make 
the special choice such that Eqs. |3.15| - |6.19 lead to fixed 
points (see Eqs. 6.36 - 6.3S ), t(Z«), n(l*), ... in Eq, 6.23j go 
to constants, independent of (and hence L±), as L± 
and hence Z», go to infinity. Thus, in this limit, we obtain 
from Eqj6~2l 



Rg(L±, L y ; t,K, . . .) — L^Rc(l, Ly/Lj_;t*, «;*, . 



(6.24) 



where t„, «;*, 



are the fixed point values of coupling con- 
stants. This result clearly agrees with the scaling forms 
for Rq, Eq.1.1 (with analogous derivation for h rms ) if we 
define Sr(xJ= Rg(l,x; t*, n*, g* n v*). 

The recursion relations Eqs.6.15-6.19 reproduce all of 



our phantom membrane results, as well as the upper 
critical embedding dimensio n df^ for self-avoidance pre- 
dicted by Flory theory, Eq.6.7, and the upper critical 



intrinsic dimension D uc — 5/2 for anomalous elasticity 
for phantom membranes. To see this, consider first the 
phantom membrane; i.e., v = 0. In this case, ft(v) = 0, 
an d to k eep t(l) fixed we see from the recursion relation 
Eq.6.15 for t(l) that we must choose 



2v + z + D-3 = 



(6.25) 



Assuming for the moment that f K {g y , g±_) — > as I — > oo, 
which, as we shall see in a moment, it does for phan- 
tom membranes for D > 3/2, we see from the recursion 
relation Eq.6.1£ for n(l) that we must choose 



2^-3z + D-l = 



(6.26) 



Solving Eqs.6.25 and 3.26 for z and v yields the phantom 
membr ane resul ts z = 1/2, v = (5 — 2D) /A, as obtained 



in Eqs.5.3J and 5.39 



To extract the upper-critical embedding dimension 
dut f° r self- avoidance from the renormalization group 
recursion relations, we construct from them a flow equa- 
tion for a dimensionless coupling constant 



vt a K b 



(6.27) 



where a and b will be chosen to eliminate the arbitrary 
rescaling exponents v and z from the recursion relation 
for v. This requirement lead to the choice 



a 

b : 



(3d-5)/8 , 
(d + l)/8, 



(6.28) 
(6.29) 



which implies 
dv 



dl 



= [6D - 1 - (5 - 2D)d] /4 - A 



" Ik 



3d - 5 



(6.30) 



Of course, an identical flow equation is obtained for v(l) 
if one instead requires that t(l) and n(l) are fixed, i.e., 
independent of I, the reby determining v and z and using 



them inside Eq. 3.19 



It is easy to see that the sign of the terms in the square 
bracket determines the relevance of the self-avoiding in- 
teraction, which becomes relevant when 



6D - 1 - (5 - 2D)d > , 



(6.31) 



i.e., for d < d^ = (6D - l)/(5 - 2D), consistent with 
the analysis of the Flory theory, Eq.6.7. 

Likewise, the renormalization group flow equations 
contain information about the upper-critical intrinsic di- 
mension for the anomalous elasticity, D uc , below which 
tubule elasticity becomes anomalous. This can be seen 
(analogously to the discussion of the releva nce of self- 
avoidance coupling v) by using Eqs.6.15-6.17 to construct 
the renormalization group flow equation for the dimen- 
sionless coupling constant 
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9y 



9v 



£3/4^5/4 ' 

chosen such that its flow 



dg. 
dl 



5 

4" 



3 f 



9y 



(6.32) 



(6.33) 



is independent of the arbitrary rescaling exponents z and 
v. Again the same recursion relation can be obtained by 
instead using the values of z and v required to kee p t(l) 
and k(1) fixed inside the flow equation for g y (l), Eq. |6.17 



It is then obvious that anharmonic elasticity becomes 
relevant for D < D uc = 5/2, where anomalous elastic- 
ity of the tubule is induced. As we will see below, in a 
phantom tubule or a tubule embedded in d > d*, this 
anomalous elasticity manifests itself only in phonon (u) 
fluctuations, i.e., softens g y , but does not renormalize the 
bending rigidity k. In physical tubules, however, which 
are self-avoiding and are embedded ind = 3<G?»~6.5, 
the elasticity is fully anomalous, both with respect to the 
phonon u fluctuations (i.e. g y vanishes as q — > 0) and the 
height h undulations (i.e. n diverges as q — > 0). 

To further analyze the renormalization of k in a self- 
avoiding membrane, it is convenient to integrate out the 
phonon field u as we did in Sec.^ for the phantom tubule, 
obtaining 



F = 



1 



2 

4- F 

~ M an 



d D - 1 x ± dy[K{d 2 y hf 
h [h] + F SA [h] 



t{dihf 



(6.34) 



where, F an h is the non-local interaction, Eq.5.11, medi- 
ated by integrated out phonons, with a kernel 



Vh(q) 



9y9±Ql_ 



9yQy 



9±Qj 



(6.35) 



and Fsa is the self-avoiding interaction. 

The long wavelength properties of the tubule phase 
will very much depend on the behavior of the denomina- 
tor in the kernel Vh at long length scales. If g y (q)q y » 
5J-(q)<7j_ ( as we saw f° r a phantom tubule) then at long 
scales Vft(q) w g±q±/q y , which behaves like ~ q y in the 
relevant limit of q± ~ qz. In this case, simple power 
counting around the Gaussian fixed point then shows 
that this elastic nonlinearity only becomes relevant for 
D < D uc = 3/2, i.e. is irrelevant for a physical D = 2- 
dimensional tubule, as we argued in Sec.|V|. 

On the other hand, if the scaling is such that gx(o)<lx 
dominates over g y (q)qy, then V/j(q) » g y , i.e. a constant 
at long length scales. Simple power-counting then shows 
that this coupling is relevant for D < D uc — 5/2 and the 
bending rigidity modulus of a D = 2-dimensional tubule 
is anomalous in this case. 

As we saw in our analysis of a phantom tubule, for 
which one is perturbing around a Gaussian fixed point 



described by q± 



q; « 



(in the long wavelength 



3/2 and k is not anomalous. We now need to extend this 
analysis to a physical tubule, i.e., to include the effects 
of self- avoidance. 

The analysis of the behavior of V/,(q) (which deter- 
mines the relevance of anharmonic elasticity) at long 
scales, around an arbitrary fixed point, is more conve- 
niently done using the language of the r enorm aliz ation 
group through the recursion relations Eqs. 5.15| and 3. IE . 
At the globally stable fixed point, in the presence of both 
the nonlinear elasticity and the self-avoiding interaction, 
we can keep t = k = 1 and g y and v fixed at fixed point 
values, by requiring 

2v + z + D-Z- f t ( v *) = , (6.36) 

2^-3z + ^-l + / K ( 5 ;, 5 l) = 0, (6.37) 

4v-3z + D-l-f g ( 9 ;) = 0, (6.38) 

2{D-l)+z- v[d - 1) - f v (v*) = . (6.39) 



In light of the above discussion, the anharmonic vertex 
for h in this renormalization group picture becomes rele- 
vant when g±(l — > oo) renormalizes to infinity, while it is 
irrelevant when g±(l — * oo) flows to zero. Thus, the rele- 
vance of Vh is decided by the sign of the r enormalization 
group flow eigenvalue of g±(l) in Eq j6.18| 



\ g± = Av - z + D - 3 



(6.40) 



which is exactly determined by the values of v and z, 
since g± suffers no graphical renormalization. 

As we have discussed in previous sections, for a phan- 
tom tubule u=(5- 2D) /A and z = 1/2. For d < d^ = 
(6D-l)/(5-2D) (= 11 fori) = 2), these values are mod- 
ified by the self-avoiding interaction, but only by order 
£ = d-d^, i.e. 



v= (5-2D)/4 + 0(e) 
z = 1/2 + 0(e) . 



(6.41) 
(6.42) 



Hence a D — 2-dimensional tubule, embedded in d di- 
mensions close to d^f = 11, \ g± = — 1/2 and g±(l) flows 
according to 



dg± 
dl 



l + 0(e)]g ± 



(6.43) 



g± is irrelevant near d 

2 /„2 2-°(V 



limit), the anharmonic nonlinearity is irrelevant for D > 



11 (for e « 1), Vfc(q) - 
g±q\/q 2 ~ q y °' e ' is irrelevant for a ph ysical D = 2- 
dimensional tubule , and, hence, f K in Eq. 3.16 vanishes 
as I — > oo. So k is unrenormalized near d — 11, for 
D = 2. That is, as we described above, the anharmonic 
elasticity is irrelevant to the bend elasticity for embedding 
dimensions near , and in this case the full model of a 
self-avoiding tubule with nonlinear elasticity reduces to 
the linear elastic truncated model introduced by usEJ and 
recently further analyzed in Ref. H. 

In this simpler (but unphysical) case, one is justified 
in ignoring the nonlinear elasticity. One is then able to 
analyze (perturbatively in e = d^c — d) t ne effects of the 
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self- avoiding interaction alone, by comp utin g the nmq =. 
tions f t (v) and f v (v) appearing in Eqs. 3.15| and 6.19| .B 



Since, as we discussed above, the "tubule gauge" sym- 
metry guarantees that in this case the self-avoiding in- 
teraction alone cann ot renormalize K, f K — 0. Thus, for 
d near df^ , Eq. |6.37 , leads to n K = and an exact expo- 
nent relation (leaving only a single independent tubule 
shape exponent): 



1 



{2v + D- 1) 



(6.44) 



which is exact for a finite range d* < d < d^c of embed- 
ding dimensions, and for phantom tubules in any em- 
bedding dimension. This result has been independently 
obtained in Ref. 0. 

However, this si mple scenario, and, in particular, the 
scaling relation Eq. |6.44 , is guaranteed to break down as d 
is reduced. The reason for this is that, as d decreases, v 
increases, and eventually becomes so large that the eigen- 
value X g± of g± changes sign and becomes positive. As 
discussed earlier, once this happens, the nonlinear ver- 
tex Eq.6.35 becomes relevant, and k acquires a divergent 
renormalization, i.e., f K ^ 0, and bend tubule elasticity 
becomes anomalous. We will now show that the criti- 
cal dimension d* below which this happens for D = 2 is 
guaranteed to be > 7/2, and hence, obviously, > 3. 

To show this, we use the exponent relation Eq. |6.44 , 
which is valid for d > d*, inside the expression for the 
eigenvalue X g± , Eq. 3.4C , obtaining 



1 



(10^ + 2L> 



(6.45) 



We then take advantage of a rigorous lower bound on v 
D- 1 



v > 



d- 1 



(6.46) 



imposed by the condition that the monomer density 

D-l-v{d-l) 



p cx L^/R^ 1 oc 



thermo dynam ic L± — > oo limit, 
side Eq.6.45 we obtain 



remain finite in the 
Using this bound in- 



A 



If D-l 

> - 10- 

~ 3 V d-l 



2D -8 



(6.47) 



from which it follows that X g± must become positive for 
d < d l *(D) with 



4 6 (2) 



4L> 



A-D 
7/2, 



(6.48) 
(6.49) 



as asserted above. 

In fact, d*(2) is probably quite a bit bigger than its 7/2 
lower bound, as two estimates of it indicate. If, for exam- 
ple, w e tak e the Flory tubule exponent v = (D+l)/ (d+1) 
in Eq. 6.45 , we obtain: 



jF 6D + 1 
df (2) = 13/2 , 



(6.50) 
(6.51) 



while if we use the e = 1 1 — ^-expansion result for v of 
Bowick and Guitter, (D — 2) 



1 



with 



we obtain 



4 + 5 2 ' 



-1.05 



5.92 



(6.52) 



(6.53) 



(6.54) 



So, based on the above estimates, we expect that in a 
D = 2-dimensional tubule, embedded in d < d* « 6, the 
fixed point of the truncated tubule model introduced by 
uscl and studied in Ref. [|, is unstable to anharmonic elas- 
ticity F an h- This means that k diverges at long length 
scales, and the scaling relation Eq.6.44 between z and v 
breaks down. Thus, for the physical embedding dimen- 
sion d = 3, the tubule bend elasticity is certainly anoma- 
lous, in the sense that k diverges, and probably quite 
strongly. We have summarized the above discussion in 
Fig.^, schematically illustrating how the renormalization 
group flow of g± , and therefore the anomalous k elastic- 
ity, change (at d*) as a function of embedding dimension 
d. 
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FIG. 6. Schematic illustration (specialized to D = 2) of 
change in relevance of g±{l) which occurs at d t . For embed- 
ding dimensions below d, (which includes the physical case of 
d = 3), g±{l) becomes relevant, leading to anomalous bend- 
ing elasticity with «(q) ~ ly 71 " j which diverges at long length 
scales. Other consequences of this qualitative and quantita- 
tive change for d < d, are discussed in the text. 
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Once d < d„, the new nontrivial relations Eqs.3.37 and 



6.38 hold, with functions f K (9y,9±) an d fg(gy) evaluated 
at the fixed point values g* and g±. 

Using the sort of renormalization group correla- 
tion fmction matching calculations described earlier, 
Eqs. 6.22 - 6.24 , it is straightforward to show that the cor- 
relation functions of the tubule, including anomalous 
elastic effects, are correctly given by the harmonic results, 
Eqs.5.6 and 5.7, except that the elastic constants g y and 
k must be replaced by wavevector dependent quantities 
that vanish and diverge, respectively as q — > 0: 



&/(q) = Qy^Sgiqy/ql) , 



with 



zrju 



fn(g* y ,gl) > 

fM ■ 



(6.55) 
(6.56) 



(6.57) 
(6.58) 



Our earlier conclusion that t he re levance of Vt is de- 
termined by the sign of X g± (Eq. |6.40| ) can be reproduced 
by simply noting that g y (c£)q^ scales like q^ +2 , and in 
the long wavelength limit is therefore subdominant to 
when 



g±ii 



2/z 
Qv 



zrj u > 2 — 2z 



(6.59) 



which, upon using Eq.6.38 and the definition of r\. 



f g (g y )/z, is identical to the condition that X g± > 0. 
The scaling functions have the asymptotic forms 



S g (x 
S K (x 



0) 
0) 



(6.60) 
(6.61) 



Combining the expressions Eqs. |6.57 and 3.58 for rj K 
and rj u with the RG fixed point conditions Eqs.6.37 and 



3.38 



shows that, at this new globally stable fixed point, 
two exact relations hold between four independent expo- 
nents z, v, r) K , and r) u (instead of a single relation Eq.6.44 
between two exponents) 



1 



-{2v + D- 1) 



3 - Tj K 
z = -J—(4 v + D-l) 

3 + ?7n 



(6.62) 
(6.63) 



That is, in contrast to the behavior for d > d*, for d < d„ 
there are two independent exponents characterizing the 
tubule phase, not one. We furthermore note that these 
exponent relations automatically contain the rotational 
symmetry Ward identity. This can be easily seen by elim- 



inating v from Eqs.6.62 and 6.63, obtaining 



2ij K + n u = 3-(D-l)/z 



(6.64) 



Ultimately, the origin of this relation is the requirement 
that graphical corrections do not change the form of the 
rotationally invariant operator (d y u + ^(d y h) 2 ). 



Just as the divergence of k is controlled by f K (gl,9±), 
the softening of g y (q) ~ q^ u is determined by the rj u = 
zf g (g*). Because f g (0) = 0, this physical g y (q) remains 
non-zero and finite as q — > 0, only if the running cou- 
pling Q v (l) in the renormalization group recursion equa- 
tion Eq. 6.17 does go to zero (because then the graphical 
piece f a ( a*) v anishes). Examining the flow equation for 



(/), Eq. 6.17 , for g y {l) to vanish, we must have 



Av - 3z + D - 1 < 



(6.65) 



However, using the lower bound on ^, Eq.6.4£ in the 
physical case of D = 2 and d — 3, we find v > 1/2. 
Hence, as long as z < 1, the condition Eq. |6.65 is not sat- 
isfied, and therefore g y (q — * 0) — ► 0, that is, r\ u > 0. We 
summarize the above discussion in Fig.^. 




4 D 



FIG. 7. Schematic of the tubule "phase" diagram in the 
embedding d vs intrinsic D dimensions. Self-avoiding interac- 
tion becomes relevant for d < d^(D) = (6D - l)/(5 - 2D), 
(= 11, for D = 2). Below the d*(Z?) curve (for which the 
lower bound is d l *(D) = {AD — l)/(4 - D)) the anharmonic 
elasticity becomes relevant, leading to anomalous elasticity 
with a divergent bending rigidity. 



We now show that the above general analysis of tubule 
anomalous elasticity in the presence of self avoidance, 
obtained using the renormalization group, can be re- 
produced via a heuristic, but beautiful physical p*gu- 
ment similar to that used by Landau and LifshitzE] to 
derive shell theory. For a tubule of diameter Rq, the 
non-zero shear g y elasticity leads to an effective -Re- 
dependent bending rigidity modulus which will be Lj_ 
and Ly-dependent if the tubule diameter depends on Lj_ 
and L y . This can be seen as follows (see Fig]|): 
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FIG. 8. Illustration of the physical mechanism for the en- 
hancement of the bending rigidity k by the shear g y elasticity. 
To bend a polymerized tubule of thickness Rg into an arc of 
radius R c requires Rg/R c fraction of bond stretching and 
therefore costs elastic shear energy, which when interpreted 
as bending energy leads to a length-scale dependent renor- 
ma lizati on of the bending rigidity k and to the Ward identity 
Eq.5.67, as described in more detail in the text. 



If we bend the tubule with some radius of curvature 
Rc Rg, simple geometry tells us that this will in- 
duce a strain e ~ d y u along the tubule axis of order 
e ~ Rq/Rc, since the outer edge of the tubule must be 
stretched by this factor, and the inner edge compressed 
by it, in order to accomplish the required bend. This 
strain induces an additional elastic energy density (i.e., 
additional to those coming from the bare k), namely 
those coming from the u elastic energy. This goes like 
g y (L y ,L±)e 2 = g y (L y , L±)(R G (L y )/ R c ) 2 . Interpreting 
this additional energy as an effective bending energy den- 
sity k v (L±, L y )/ R 2 , leads to the effective bending mod- 
ulus K y (L± 7 L y ), 



K y (L±,L y ) ~ g y (Lj 



> L v )Rg(Lj 



(6.66) 



Inserting the scaling forms n y (L±,L y ) = L^ K S K (L y /Lj_), 
g y {L ± ,L y ) = L-^S g (L v /L z ± ) and R G (L ± ,L y ) = 
Lj_Sii(Ly/Lj_) into above expression, we obtain a rela- 
tion between the scaling exponents 



2v = z(j] K + rj u ) 



(6.67) 



which is exactl y th e expone nt re lation one obtains by 
subtracti ng E q. |6.37 from Eq.6.3£, and using the expres- 
sions Eq.6.57 and |6.5§| for r/ K and rj u , all of which were 
obtained using renormalization group arguments. 

Since t he a bove phys ical shell argument is very gen- 
eral, Eqs. 5.66 and 6.67 hold independent of the mecha- 
nism that generates anomalous elasticity. For th e case 
of the phantom membrane (for D > 3/2) Eq. |6.66| reveals 
that k is not anomalous because the softening of the shear 
modulus gi/(q) by thermal fluctuations precisely compen- 
sates the bending rigidity pr oduce d by the finite diameter 
Rg of the tubule. Equation 6.67 then correctly predicts 



for the phantom tubule that r\ u — 2v/z, which is con- 
sistent with the phantom tubule results rj u = 5 — 2D, 
v = (5 — 2D) J A, and z = 1/2. Furthermore, because the 
anharmonic elasticity Vh(q) is irrelevant for d> d„, 

r)u = 2u/z , (6.68) 

is valid, even in a self-avoiding tubule embedded in these 
high dimensions. 

We note, finally, that all of the exponents must show 
a jump discontinuity at d*, as shown in Fig.^. There- 
fore, unfortunately, an extrapolation from e = 11 — gL. 
expansion in a truncated model with linear elasticityu 
down to the physical dimension of d — 3 (which is below 
<i*) gives little information about the properties of a real 
tubule. 



exponents 




FIG. 9. Schematic graph of the shape exponent v and 
anomalous bend exponent g K (for D = 2). Note the jump 
discontinuity as a function of embedding dimension d, occur- 
ring at d = d* ~ 6. 

The computations for a physical tubule must be per- 
formed for d < d*, where both the self-avoidance and the 
anharmonic nonlinearities are relevant and must be han- 
dled simultaneously. As we discussed above, for d < <i„, 
the eigenvalue X g± > 0, leading to the flow of g±(l) to 
infinity, which in turn leads to Vft,(q) = g y . Physically 
this regime of g± — + oo corresponds to freezing out the 
phonons u, i.e. setting u — in the free energy F[h, u] 
in Eq. 6.13| . This is consistent with our finding that for 
d < d„, in the effective free energy F[h] (with phonons in- 
tegrated out), Eq. 3.34, the kernel Vh = g y - The resulting 
effective free energy functional for a physical self-avoiding 
tubule is 



x±dy 



K(d 2 y hf +t{8^hf +- A 9y{d y h) 



+ v IdydP-^dP-^SWfiixx^-K&^y)) , 

(6.69) 
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Unfortunately, no controlled perturbative study is possi- 
ble for d < d*, since one must perturb in g y around a 
nontrivial, strong coupling fixed described by v* — 0(1) 
and g* = 0. Furthermore, as we will show below, at 
this fixed point there is no upper critical dimension for 
g y , i.e. anharmonic nonlinearities are always relevant for 
d = 3 < d*, for any D. This strongly contrasts with the 
Gaussian fixed point (describing phantom membranes) 
at which the anharmonic nonlinearity is only relevant for 
D < D uc = 5/2. 

In what follows, we will illustrate how one might at- 
tempt to actually calculate the exponents v, z, t] k , and 
T] u , for d < d*, and enumerate the (many) technical diffi- 
culties that prevent us from doing so, and conclude with 
a cautionary list of several unsuccessful uncontrolled ap- 
proximations that we have tried. 

In principle, all we need to do is calculate t he h 
(i = t, v, g, k) functions in the recursion relation Eqs.6.15- 



6.19 , which represent the perturbative ("graphical") cor- 
rections to the associated coupling constants. Once these 
/-functions are known they give 4 equations (Eqs.6.36- 
6.39) that uniquely determining the 4 unknowns tubule 
shape exponents, v, z, v* , and g*, as well as the the 
flow of g±(l), and therefore completely characterize the 
long wavelength properties of self-avoiding anharmonic 
tubules. 

Our goal then is to calculate /t(f), fv{v), fn{9y), and 
f g {g y )- The functions f g {g y ) and f K (gy) are determined 
by the diagrammatic corrections to g y and k, with the 
corresponding Feynman diagrams displayed in Fig.[ic|. 
The results to leading order in g y , are 



fnidy) = C K g 2 y , 
fg\9y) — C g g y , 



(6.70) 
(6.71) 



where C' K and C g are d and .D-dependent constants, 
whose calculation proves to be the sticking point, as we 
will describe below. 




(a) 




FIG. 10. Feynman graphs that renormalize: (a) the anhar- 
monic elasticity g H ,and (b) the bending rigidity k. 



Of course, once d is below no matter how close it is 
to d* , the fixed point that controls the elastic properties 
of the tubule phase is not perturbative in g y . That is, we 
do not expect g y to be 0{d* — d), but, rather, O(l), even 
for d* — d << 1. Furthermore, of course, since <i* ~ 6, 
d* — d is not small in the physical case d = 3 anyway. 
For both of these reasons, truncating the calculations of 
f K and f g at the leading order in g yi as we have done in 
Eqs.|6~70|and|6Jl, is an uncontrolled, and far from trust- 
worthy approximation. However, we know of no other 
analytical approach. Furthermore, as we shall see, even 
this uncontrolled analytic approach proves intractable: a 
reliable calculation of the values of the constants C K and 
C g has eluded us. 

To complete the characterization of the fixed point we 
can proceed in two ways. The most direct way is to 
simply pcrturbatively evaluate the functions ft(v), fv(v). 
Luckily (foij-|Us) this has recently been done by Bowick 
and Guittem in a truncated harmonic tiibule model (pre- 
viously introduced and studied by usEl) near d = d^. 
Although, for the reasons that we discussed above, these 
calculations are not rigorously applicable to a physical 
tubule in d = 3 < d* (where anharmonic elasticity is cer- 
tainly important), for lack of being able to do any bettes 
we extrapolate these functions, computed near d = 11 B 
down to d = 3 



ft(v) - C t v , 
f v (v) = C v v , 



(6.72) 
(6.73) 



Now using Eqs. |6.70| - |6.73 in Eqs. 3.36 - 5.39| we obtain 
four equations for four unknowns (z, v, g*, and v*), ex- 
pressed in terms constants C K , C g , Ct and C v , (special- 
ized here to D = 2). 



2u + z-l- C t v* = 
2v-3z+l + C K gf = 

4u-3z + i- c g9 ; = 

2 + z-v[d-l)- C v v* = 



(6.74) 
(6.75) 
(6.76) 
(6.77) 



where the constants Ct and C v (computed in the trun- 
cated tubule model near d = d^ for D = 2) are given 
byEl 



C t 



C v 



8ir 2 ' 
0.068 



(6.78) 
(6.79) 



These equations can be uniquely solved for v, z, g y , and 
v*. In terms of C K and C g , in D = 2 and d = 3 we obtain 
for v and z 



_ 1 1 Cg 

Z ~ 3 + 6C K + 6C K 



(6.80) 
(6.81) 
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from which rj K and r\ u ca n also be deter mined using the 
solution for g* inside Eqs. 6.57 and 6.5§| 



3C 



Vk = 



1 + C q + 2C K ' 



3-3C 



1 + C„ + 2C K 



(6.82) 
(6.83) 



Another approach to estimating the tubule shape ex- 
ponents is to rely on the usual accuracy of the Flory the- 
ory (in treating the effects of self-avoidance), instead of 
the extrapolation of functions f v (v) and ft{v) down from 
e-expansion. Although it is usually not stated this way, 
in the language of renormalization group, Flory theory 
amounts to assuming that the graphical corrections to t 
and to v are the same, i.e. /„(«*) = ft(v*). Using this in 



Eqs.6.36 and 6.39, we obtain the Flory result for v 



D + l 



VF 



= - . for d = 3 , D 

4 



(6.84) 
(6.85) 



consistent with our earlier analysis in subsection VI A . 
Note that, if f v (v) — ft(v) for all v, this result would be 
exact independent of the jump in the other exponents z, 
rj K , and rj u at d*. That is, it would apply even below d*, 
and v would not jump, or be in any way non-analytic, at 
d*. 

Now, of course, we know from the explicit leading order 
calculation in Ref. || that f v (v) does not = ft(v) exactly. 
However, we do know from that calculation that they are 
quite close, at least to leading order, as illustrated by the 
good agreement between Flory theory and the extrapo- 
lated e-expansion. If this persists down to d = 3, and 
to large v, and our experience with polymers suggests 
that it will, then v may be quite accurately predicted by 
Flory theory, despite the complications associated with 
the onset of anomalous bend elasti city a t c?* . 

U sing the Flory value for v (Eq. |6.84 ) inside Eqs. 5.75 
and 6.76, together wi th th e di agram matic corrections to 
K and q y given in Eqs.|/7(] and we obtain two equa- 
tions (specialized to D = 2) 



6/(d+l)-3z + l + C* K . 9l * =0, 
12/(d+l)-3z + l-C 9 . g ; = 0, 



which gives for d = 3 



4 
3 



r<2 

6C K 



C fl (C» + 6C H )V2 



6C, 



(6.86) 
(6.87) 



(6.88) 



Now, at least in this uncontrolled approximation of 
truncated perturbation theory at one loop order, it seems 
that we are left with the straightforward task of calcu- 
lating the constants C K and C g . Alas, things are not so 
simple, for reasons that are undoubtedly connected with 
the fact that d* is not perturbatively close to , which 



is the only dimension about which one can do a gen- 
uinely controlled approximationBH, and the much more 
surprising fact that, even though ej = 5/2 — D is only 1/2 
(for D = 2), this e/-expansion in intrinsic dimension, as 
we will show, is demonstrably extremely unreliable, giv- 
ing qualitatively different answers, such as a reduction, 
rather than an increase of k due to fluctuations. 

Our unsuccessful (but heroic) attempts to calculate C g 
and C K were as follows: 

(I) 

Calculate them in an ej = 5/2 — D-expansion for a 
phantom membrane, then use these same constants C g 
and C K for the real, self-avoiding membrane. This ap- 
proach obviously makes many errors, since, by the time 
we get down to (5/2), the correlation functions of the 
true, self-avoiding membrane are already quite different 
from those of the phantom membrane, due to the effects 
of self-avoidance. Furthermore, these effects are particu- 
larly pronounced for intrinsic dimensions D = 5/2, since 
d$£(5/2) = ao, as illustrated in Fig.0. 

Nonetheless, since no other analytical calculation is 
available (and we are persistent young lads), we at- 
tempted this e/ = 5/2 — D-expansion. However, the 
results made no physical sense: we found a negative r] K , 
i.e., a downward renormalization of k. The detailed cal- 
culations are virtually identical to those for the renor- 
malization of k at the tubule-to -cru mpled phase transi- 
tion, which are described in Sec. VII. We note here sim- 
ply that the origin of this negative contribution to k is 
a negative region of the real-space correlation function 



G(x ± ,y) = Cix^/y^^Yix^/y 2 ), as given by Eg .[731 



The integrand x 5 ' 4 Y 3 (x) in the x- integral of Eq. 7.37 has 
a negative region which, though narrow, actually over- 
whelms the positive contribution to r\ K from the much 
longer, but smaller, jtail, as we have verified by direct 
numerical integratior£3. 

This negative region is purely an artifact of calculating 
in a fractional intrinsic dimension D = 5/2. In D = 2 for 
a phantom membrane, where there is no relevant anoma- 
lous elasticity for h, and hence we can calculate ft.- ft. cor- 
relation functions exactly, we find the analog of Eq. 7.31 



G(x,y) = 



dq x dq y 



Jq^x+iqyy 2 



(2tt)2 
1 



4(tt|x|) 1 /2 



q. 

-y 2 /(i\x 



4 



(6.89) 
(6.90) 



which, unlike the analogous correlation function in D = 
5/2, Eq. 7.31 , is positive definite. Thus, the anomalous 
contribution to k in D = 2 will also be positive, as we ex- 
pect on physical g round s (i.e., the shell theory argument 
summarized in Eq. |6.66D , while the 5/2 — D-expansion is 
qualitatively wrong in predicting a negative renormaliza- 
tion of k. Clearly, it cannot be trusted quantitatively 
either, and is, in fact, totally useless. 
(II) 
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Direct, uncontrolled RG in D = 2. Now, we at least get 
qualitatively correct upward renormalization of k. How- 
ever, here we have a different problem, that appears in 
any perturbative calculation away from an upper critical 
dimension (and is usually "swept under the rug"): even 
though D = 2 would not, a priori, appear to be far below 
D = 5/2, it is, in the sense that graphs that only diverge 
logarithmically in D = 5/2 diverge extremely strongly in 
D = 2. In particular, f ollow ing very closely the manipu- 
lations that lead to Eq. 7.35] , we find a contribution to k 
of the form 



5n 



Cl 



ydy 



dx- 



-3/(4x) 



(6.91) 



where c\ is a well-determined constant that we could cal- 
culate, and C2 is an arbitrary constant which depends on 
precisely how the infrared divergence of the above inte- 
gral is cutoff with q y . This a rbitrary constant is the prob- 
lem: if the integral Eq,6.91 had diverged logarithmically, 
the precise value of the constant C2 would be unimpor- 
tant (it would just lead to a finit e additive constant). 
But, since the integral in Eq. 6.91 diverges so strongly 
(like (c2/q y ) 2 ) in D = 2, it is extremely sensitive to the 
precise value of C2, which we have no clue as how to 
choose. Thus, we have no ability to predict r\ K at all by 
this approach. 

This strong divergence indicates that in this sense 
D = 2 is quite far from D = 5/2, and any kind of per- 
turbative approach, even to simply calculating one loop 
constants like C g and C K) is doomed. 



C. Gaussian variational theory of self-avoiding 
tubules 

Here we study the effects of self-avoidance within the 
tubule phase using the Gaussian variational method, 
which was previously applied to the st udy .of self- 
avoidancp-in crumpled isotropic membranesELrEa and in 
polymers£-l. It is important to emphasize that both Flory 
theory and the Gaussian variational method are uncon- 
trolled approximations in that there is no way to system- 
atically estimate and reduce the error. 

We begin with the effective Hamiltonian that describes 
the long wavelength behavior of the tubule for d < d*. 



H = 



1 



jD-l 



x±dy 



K (d 2 h) 2 +t(d^h) 2 + -g y (d y h) 4 



+ v Jdyd D - 1 x^d D - 1 x' x 5^- x \h{^y)-h{^y)) , 

(6.92) 

where, in contrast to other sections we use the notation 
H to distinguish long-wavelength effective Hamiltonian 
(the free energy functional) from the actual free energy 
F . Computation of correlation functions in the presence 
of the self-avoiding nonlinearity cannot be done exactly. 



However, we can replace the Hamiltonian H, Eq.6.92 



by a variational Hamiltonian H v , quadratic in the fields 
h(x±,y), which allows exact calculations of any corre- 
lation function. Following the standard variational pro- 
cedure, we then pick the "best" form of this variational 
Hamiltonian, where by "best" we mean thakit minimizes 
an upper bound on the true free energy F'SB 



F<F=(H-H V ) V +F V 



(6.93) 



We take our variation ansatz Hamiltonian to be 



H„ 



dkyd - 1 ^ 
(2tt) d 



G V {k ± ,ky)\h(k X ,ky)\ 2 , (6.94) 



where G v (k±, k y ) is the variational kernel to be opti- 
mized over. Note that because of anisotropy intrinsic to 
the tubule, G v (k±, k y ) is not rotationally invariant as it 
is for the analogous analysis of the crumpled phase. 



We now compute the right-hand-side of Eq.6.93 and 
minimize it over G„(kj_, k y ). 



(H - H„ 



A 



G V (k ± ,ky))(\h(k ± ,k y )\ 2 ) v 



(6.95) 



where A = L»,L , 1 is the "area" of the membrane 



and we defined J = Jdyd D x x±d 



_L' 



and / 



/ dkyd^kx /(2tt) d . 

The above averages are easily evaluated with 
(\h(k ± ,k y )\ 2 ) v = (d - 1)/G„(k) and (S) v = 
(^-D (h(x ± ,y) - h(x' ± ,y))) v given by 

IXS - I f ^ % A r i ?-(Hx±,y)-H*', ,y))\ 



d d - l q 
{2-Kf- 1 

1 



-g 2 K(|xx-x'jJ) 



(27T) 



d-1 



K(\x ± -x> ± \) 



(d-l)/2 



(6.96) 



where, 



K(\x ± \) 



1 



-(|fc(x 1> Q)-/i(Ox ) 0)| s 



2{d- 1) 

1 — cos(kj_ • xj_)J 
/ ; GJk) 



(6.97) 



and we have used in Eq.6.96 the Fourier representation 
of the d — 1-dimensional delta- function. 

Putting all th is together we obtain for the right-hand- 
side of Eq]6~93l 
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- 1 



9v(d+l) 
4(d-l) 



k 2 



{2n) d - 1 {d-l) 

+ [ io g (a(k)) , 

Jk 



l k G v (k) 

v (d-l)/2 



that the right hand side of Eq.6.102 is proportional to 
xT with 



7 = 1 - D + 3a , 

_ 7(1 - D) + 3v(d + 1) 



(6.103) 



(6.98) 



which when minimized with respect to G v (k), 
SF/6G v (k) — gives an integral equation 



G v (k) = nki + tk]_ - 



2v 



(4 7r )(d-l)/i 



iD-1 



K ( x± yd+i)/2 



(6.99) 



The only effect of the anharmonic elasticity term g y is 
to generate an upward renormalization of the effective 
tension along the y-axis 



Sty = 



9y(d+l) 
2(d-l) J k G v (k) 



(6.100) 



Since we must choose the renormalized tension along the 
extended tubule axis (y) to be exactly zero in order to 
treat the free tubule, all of the anharmonic elastic effects 
disappear in this Gaussian variational approximation. 
That is, to correctly model a tubule with free bound- 
aries, we should have started with an elastic Hamiltonian 
with a bare, negative tension piece that exactly c ancell ed 
thermally generated positive contribution in E q.|6.100 , 



The simultaneous integral equations Eq. 3.97 and 
Eq. |6.99 determine G v (k) and K(x±). At long length 
scales they ar e sol ved by K(x±) ~ x", where from the 
definition Eq.6.97, we see that K(x± = L±) is propor- 
tional to the square of the radius of gyration or the tubule 
thickness that we are after, and hence the v that solves 
these coupled non-linear integral equations will be the 
Gaussian variational prediction for the radius of gyra- 
tion exponent as well. We substitute this scaling ansatz 



into Eq. 6.9£ for G v (k), and find that while for d > d^c 
the self-avoidance is irrelevant and v = (5 — 2D) /A (as 



found in Sec. VI A), for d < df^ these int egral equations 
can only be solved if the tk\ term in Eq. 6.99] is exactly 
cancelled by a part coming from the integral in the last 
term and the resulting propagator takes the form 



G„(k) = Kk\ 



D+l 



(6.101) 



where v cx v is an effective self-avoiding inter action pa- 
rameter. Substituting this form into Eq.6.97 for G„(k), 
and requiring self-consistency with our original ansatz 
K(x ± ) 



x 2 ^ gives 



d D 1 q±dq y (l - cos(kj_ • xj_)) 



ik ( d+ i),-D + i 



(6.102) 



Making the change of variables q± = 
q y = q y /\x±\ a , with a — (u(d + 1) — D 



cLl/|x_l| and 
f l)/4 reveals 



To satisfy the self-consistent condition Eq. |S.102 



must be equal to 2v. The resulting simple linear equation 
for v has a solution (for d < d^(D)): 



x± (1 — cos(kj_ • Xj_)) 



7D-7 

3d -5 ' 



which for the physical case of D — 2 gives 
7 



3d -5 
1 

4 



for d < 11 , 



for d > 11 



(6.104) 

(6.105) 
(6.106) 



We observe that v(d = 4) = 1, and therefore (ac- 
cording to the Gaussian variational approximation) the 
tubule is no longer crumpled along the _L-direction. This 
suggests that the tubule phase is unstable to the flat 
phase in embedding dimensions d < 4 (which unfortu- 
nately includes the physical case of d = 3). However, as 
discussed in the Introduction, the Gaussian variational 
method is an uncontrolled approximation. It probably 
does give the correct trends of, e.g., exponents with di- 
mensionality d. However, the variational approach is very 
close, in spirit and technically, to the large d expansion 
methods, and therefore intrinsically unable to obtain the 
small d dependence correctly. It is therefore difficult to 
place any faith in the actual values of exponents, partic- 
ularly when the value of v at small d actually determines 
whether the tubule phase survives or not. 

We believe that this Gaussian variation theory is in- 
correct in predicting that the tubule phase does not exist 
in the presence of self-avoidance in d = 3, and-, reiterate 
our earlier observation that both Flory theoryu and the 
e = 11 — d-expansiorJ3 predict that the tubule phase sur- 
vives self-avoidance. Since both these latter approaches 
agree quite closely with each other, and since, further- 
more, the e-expansion is the only controlled approxima- 
tion, we are far more inclined to trust them than the 
uncontrolled Gaussian approximation, which agrees with 
neither. 

The final determination of whether or not the tubule 
phase survives self-avoidance will, or course, rest upon 
simulations and experiments, both of which we hope our 
analytic work stimulates. 



VII. FLUCTUATION EFFECTS AT 

CRUMPLED-TO-TUBULE AND 
TUBULE-TO-FLAT TRANSITIONS 

The transition from the crumpled-to-flat rai^ase in 
isotropic membranes has been previously studiedcJ and is 
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predicted to be driven first order by fluctuations for em- 
bedding dimensions d < d c = 219. As can be seen from 
Fig.|l| this direct transition is very special for anisotropic 
membranes. It is easy to see that any path finely tuned to 
pass through the tetracritical point will undergo a direct 
crumplcd-to-flat transition identical to that of isotropic 
membranes, discussed in Ref. ^4|. 

Here we focus on the new transitions crumpled-to- 
tubule and tubule-to-flat, which are generic for mem- 
branes with any amount of anisotropy. As we discussed 
at the end of Sec. Ill, there are two possible mean field 



phase diagram topologies depending on the values of mi- 
croscopic elastic moduli of the membrane. However, for 
the crumpled-to-tubule transition there is no difference. 
In this section we first study the crumpled-to-tubule tran- 
sition for a phantom membrane using a detailed renor- 
malization group analyses. We then study both the 
crumpled-to-tubule and tubule-to-flat transition using 
scaling theory, incorporating the effects of both the an- 
harmonic elasticity and self-avoidance. We postpone the 
more technically challenging renormalizatiop-group anal- 
ysis of the phantom tubule-to-flat transitioned and renor- 
malization group analysis of crumpled-to-tubule and 
tubule-to-flat transitions for self-avoiding membranesES 
for future publications. 



A. Renormalization group analysis of 
crumpled-to-tubule transition 

We s tart out with the general free energy defined in 
Eq. [2.l| , for now ignoring the self-avoiding interaction. 
Without loss of generality we will study the transition 
from the crumpled to the y-tubule phase. As discussed 
above, in mean-field theory, this transition occurs when 
t y — > from above, while ij_ remains finite and pos- 
itive. Hence, simple power-counting on the quadratic 
part of the free energy leads to anisotropic scaling at the 
transition with q± oc q^. Therefore, the only relevant 
terms quadratic in r near the transition are: the bending 
rigidity along the y-direction (k v (c^r) ), and the surface 
tension terms along the y and _L-directions t y (d y r) 2 and 

t± ipa^) 2 1 respectively. The corresponding noninteract- 
ing propagator at the transition is 



The critical properties of the crumpled-to-tubule transi- 
tion can be obtained by applying scaling theory and the 
renormalization group to this free energy exactly as we 
did earlier in treating fluctuations in the tubule phase 
itself. In this case, "lengths" means intrinsic coordinates 
x = (xj_ , y) , and the "fields" are the extrinsic positions 
r(x). Because of the strong scaling anisotropy of the 
quadratic pieces of the free energy, we rescale xj_ and y 
anisotropically: 



x^ = x 

y 



y'e* 



and rescale the "fields" according to 



f(x) 



(x') 



Under this transformation 



D (D-l-3z+2 x )l 



Ky(l) = 

t ± (l) = t ± eS D - 3+z + 2 ^ 1 , 



(7.3) 
(7.4) 



(7.5) 



(7.6) 
(7.7) 



Requiring that both k v and t± remain fixed under this 
rescaling (zeroth order RG transformation) fixes the 
"anisotropy" exponent z and the "roughness" exponent 
X (which is the analog of v for the tubule phase): 



1 

Z = — , 

2 ' 

X = (5/2-D)/2 



(7.8) 
(7.9) 



Alt hou gh this choice keeps the quadratic (in r) part of F 
Eq.7.2) unchanged, it does change the quartic piece: 



Uyy {I) 



; (D-l-32+4 X )i 
»(5/2-D)l 



(7.10) 
(7.11) 



where in the second equality we have used the results 
Eqs.7.8 and 7.9 for z and \. We see that for D < 5/2, 



'mi 



(n(q)^(-q)) 



txq± + t v q% + Kyq^ 



The anisotropic scaling dictated by this noninteracting 
propagator at the transition (t y = 0) leads to significant 
simplification of the interaction term in the free energy. 
Keeping only the dominant nonlinearity we obtain 



F[r(x) 



x±dy 



grows upon rescaling. Physically, this means that its 
effects become more important at longer length scales. 
At sufficiently long length scales, it completely invali- 
dates the harmonic elastic theory and the naive pertur- 
bation theory in the nonlinearity u yy around it, even 
for arbitrarily small coupling u yy . Simple additional 
\ C V y j ( -1) anisotropic rescaling of xj^ = cac' ± and y — (3y' , with 
(3 = (tj_/Ky) 1 ' 2 a 2 , which rescales n y and tj_ to 1, reveals 
that the effective coupling consta nt of the nonlinearity is 
Uyy/ny. This, together with Eq j7.1l| , predicts that the 
characteristic length scale L™' beyond which the dimen- 
sionless coupling constant becomes of order 1 and the 
harmonic elastic theory and perturbation theory (around 
it) break down is 



+ t y (dy^ 2 + ^(dyr-d y rr 



(7.2) 



L i = 



1/(5/2-13) 



(7.12) 
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To analyze the new behavior that prevails on even 
longer length scales requires a full-blown renormalization 
group analysis. .— . 

Such an analysis^ will lead to corrections to the sim- 
ple rescaling of «„, t±, and t y , due to the non-linearities 
(in this discussed above). These corrections 

can be absorbed into "anomalous" exponents r\ K , rjt, and 
66, defined by the large renormalization group "time" 
(I — ► oo) limits of K y (l), t±(l), and t y (l), respectively: 



1 



JD-\- l iz+zri K +2 X )l 



K y (l) = 

, (1\- f (D-l-z-Se+2 x )l _ , 

y v / — y — y 



Xtl 



(7.13) 
(7.14) 
(7.15) 



The exponent At defined above is the thermal eigen- 
value of the reduced temperature (surface tension along 
y-direction) which is an inverse of the correlation length 
exponent along the _L-direction (see below). Requiring 
that K y and t± remain invariant under the renormal- 
ization group transformation determines the values the 
anisotropy exponent z and the field rescaling exponent 
X, 



2-Vt 
4-77 K ' 

10-4£> + 7 ?B (D-3 + ^)-3^ 
8-2 Vk 



(7.16) 



which, as quoted above in Eqs.7.8 and |7.9| , reduce to 
z = 1/2 and x = (5/2 — D)/2, for r] K = r\t = 0, as is valid 
at zero order in perturbation theory in u yy . 

Once the values of rj t , fj K and \ a t the critical point 
are determined, the renormalization group gives a rela- 
tion between correlation functions at or near criticality 
(small t y ) and at small wavectors (functions that are dif- 
ficult to compute, because direct perturbation theory is 
divergent) to the same correlation functions away from 
criticality and at large wavevectors (functions that can 
be accurately computed using perturbation theory). For 
example the behavior of the correlation lengths near the 
transition can be deduced in this way: 



U(ty) = e l U(t v e x * 1 ) , 



(7.17) 
(7.18) 



(7.20) 



At 2(2 -m- 266) - Vk (2 - Vt - 66) ' 
v y — zv± . (7-21) 

We now compute the anomalous exponents to lowest 
non-zero order in e, where e = 5/2 — D. As usual in 
the e-expansion, the order at which a given graphical 
correction enters the perturbation theory is equal to the 
number of loops in the associated Feynman graph. We 
split the field r(x) into short and long wavevector parts 
r (x) = r < (x) + r> (x) and integrate over the fast fields 
r>(x). Diagrammatically this leads to one-loop correc- 



tions to Uyy and t y . 



There are no corrections to K y to 



first order in e, i.e., r\ K = 0(e ). Furthermore, since 
the interaction u yy always carries a factor of q y with ev- 
ery field r, the t± tension remains unrcnormalized, and 
rjt = to all orders, implying z = 1/2 + 0(e 2 ) and 
X =(5/2-£)/2 + 0(6 2 ). 




(a) 




(b) 




(c) 



FIG. 11. Feynman graphs that renormalize: (a) the non- 
linearity Uyy, (b) the tension t y and (c) the bending rigidity 



where in the above we assumed that a critical fixed point 
exists and all other coupling constants have well-defined 
values at the fixed point. Using the above equations for 
t„e Xtl « 1, we obtain, 



U(t y )~at-^ , 
£ y {t y ) ~ a t y v , 

where a « £(1) is the microscopic cutoff and, 



(7.19) 



The first two diagrams in Fig.|TT|, followed by the rescal- 
ing introduced above (necessary to restore the original 
uv cutoff), lead to the one-loop recursion relations for 

u = (K3/2/ ' V2)u yy / (Ky^ 4 t 3 / 4 ) and t y , respectively, 



— = eu- {d + 8)u 2 , 
Bt 

-^ = (l-(d+2)u)ty, 



(7.22) 
(7.23) 
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where e = 5/2 — D and -K3/2 is the surface of area of a 
3/2-dimensional sphere divided by (2tt) 3 ^ 2 . As usual, in 
the above, we also redefined t y to be the reduced tem- 
perature, measured from its true value at the transition 
(which in mean-field theory starts out at 0, but is shifted 
to a negative value by fluctuations). Note that, in con- 
trast to the familiar e = 4 — D- expansion for critical 
phenomena, for which e = 1 in the physical case D = 3, 
here we have e = 1/2 in the physical case D = 2. Hence, 
our one-loop expansion should be quantitatively more ac- 
curate by a factor of e~ 2 = 4, than the 4 — D-expansion 
at the same order. Thus, we expect our one-loop values 
for v± and v y to be ac curate to ±0.02. 

Examining Eq.7.23, we observe that for D < D uc = 
5/2, (i.e., for e > 0) the Gaussian fixed point is unsta- 
ble and the critical properties of the crumpled-to-tubule 
transition are characterized by a nontrivial fixed point 
with a fixed point value it* of u given by: 

(7.24) 



d + 8 

Note that, in contrast to the treatment of crumpled-to- 
flat transition in isotropic membranes 2 .!, where the criti- 
cal point was only stable for an unphysically large value of 
the embedding dimension d > 219, the critical point char- 
acterizing the crumpled-to-tubule transition found here 
is stable for all d. 



Equation 7.22 c an b e easily integrated once the fixed 
point value u*, Eq.7.24 , is inserted for u; comparison with 
the general Eq.7.16 then gives At: 



A, = 1 



d + 2 
d + 8 



(7.25) 



which, upon using Eqs. 7.20] and 7.21, gives for a physical 
membrane (D — 2, d = 3) 



1.227 , 



v y fa 0.614 



(7.26) 
(7.27) 



The i] K exponent to 0(e 2 ) is determined by the dia- 
gram in FigJTT]c. Evaluating this diagram in real space 
and then Fourier transforming, we find that this con- 
tributes to the free energy 

SF = -16u 2 (d + 2) f <j2|r(q)| 2 r(q) , (7.28) 

J Q 



where 



T(q) = J d 3 / 2 x^dy e 4qx G 3 (x) , 



(7.29) 



with, in turn, 
G(x ± ,y) : 



d 3 / 2 q ± dq y e 1 ^-^ q 2 



(2tt) 5 /2 

where we have rescaled lengths so that n y = t±_ = 1 



(7.30) 



After a contour integral over q y , and an angular inte- 
gral /^^(sinfl)^- 2 '^^^, we obtain, 

G(xx,y) = 2-/V- 3 /V 2 (^) V V(^) , (7.31) 



where we have defined 



Y(x) 



duu 1/4 J_ 1/4 (xu) e-\^cos(v/u72 + 7T/4) . 

(7.32) 



Now going back to Eq. [7.28| we observe that the q = 



piece of T(q) contributes to the gj|r(q)| part of F, ad- 
ditively renormalizing t y which corresponds to the usual 
inconsequential T c (critical tension) shift. The order q 2 
piece of T(q) renormalizes K y . We define 



r( q ) = r(o) 



2<% B {qv 



(7.33) 



where 
B(q y ) 



dyd 3 / 2 x ± y 2 G 3 (x ± ,y) . (7.34) 



Note that the infra-red cutoff on the integral over y is 
q^ 1 - This integral diverges logarithmically as q y — * 0. 
We can identify the coefficient of the logarithm with r\ K 
in the expression n y (q y ) oc q y v "- 

To extract this logarithmic divergence, we make a 
change of variables in the integral |xj_| = xy 2 and find 



B(q v ) 



1 



2 13/4 7r 3/2 r (3 / / 4 ) 



f 9v d JL f dxx b ' i Y 3 {x) , 
A-i V Jo 

(7.35) 



where we have used the fact that the surface area of a 
3/2-dimensional sphere is 27r 3 / 2 /r(3/4), and taken into 
account the factor of 2 coming from the fact that the orig- 
inal integral over y extends over both y > and y < 0. 

Putting all of the above together and evaluating the 
coefficient of the \og(q y ) at the fixed point value of g* 
from Eq.7.24, we obtain 



C(2)(rf + 2) 
8(d + 8) 2 ' 



(7.36) 



where 



C(2) = 2 23 / 4 r(3/4) / dxx 5 / 4 Y 3 (x) , (7.37) 



The value of C(2) has been calculated numerically^ to 
be C(2) w -1.166±0.001. Using this value, e = 1/2, and 
d = 3 in Eq. 7.36| , we find that rj K is very small, 



r? K ( J D = 2,rf = 3)w -0.0015 



(7.38) 
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As noted earlier in our discussion of the tubule phase 
itself, we do not trust this negative value of 7] K , but, 
rather, believe it to be an artifact of the peculiar negative 
regime that appears in the correlation function G(xj_,y) 
in D — 5/2. We expect r\ K to be positive, but still quite 
small, at the 'phantom tubule-to-crumpled transition. 

Given the smallness of r\ K and e, and the vanishing of 
rjt, the exponents computed here to first order in e are 
expected to be very accurate. 



which are always crumpled, nonetheless assume, on aver- 
age, non-spherical shapescS, as can be seen, for example, 
by looking at the ratio of the average maximum and min- 
imum eig enval ues of the moment of inertia tensor. Our 
result Eq. 7.39 for t y > is only a little more surprising, 
since it predicts an aspect ratio R y /Rc that actually di- 
verges as T — ► T^., and membrane begins to extend into 
a tubule configuration. 

The exponents 1^_'J_ defined in above equation obey 
the scaling laws 



Scaling theory of crumpled-to-tubule and 
tubule-to-flat transitions 



G,y 

7+ 



,G,y 



(7.40) 



We will now incorporate the effects of self- avoidance 
on these transitions. We have not yet done the full 
renormalization group analysis of this problem (which 
must include bnih the elastic and self-avoiding interaction 
nonlinearities)t3, and limit ourselves here to discussing 
scaling theory and the Flory approximation. 

Near the crumpled-to-tubule transition, for square 
membranes of internal size L, we make the following gen- 
eral scaling ansatz for the extensions R y and Rg of the 
membrane along and orthogonal to the tubule axis, re- 
spectively: 



Rg, 



oc 



\ty\ 



ty >0,L» £ ct 

l « ict 

, ty <0,L» get 



(7.39) 
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G,y 



G,y 



G,y 



(7.41) 



As always, these scaling laws follow from requiring that 
the generalized scaling form matches on to known results 
in the appropriate limits. 

From Flory theor y we c an derive the values of the crit- 
ical exponents in Eq.7.39, as we have already derived the 
exponents v t and v c characterizing the tubule and crum- 
pled phases, simply by being more careful about tem- 
perature dependent factors in that derivation. Again, we 
start by estimating the total self-avoidance energy Eq |2.1| 
in the tubule phase (i.e., t y < 0) as Esa ~ Vp 2 . Now, 
however, we very carefully write the volume V in the em- 
bedding space occupied by the tubule as V sa R ( ^~ 1 R y . 
Writing 

Ry = (yLy , (7.42) 



where subscripts t, c and ct refer to tubule, crum- 
pled and tubule-to-crumpled transition, respectively, and 
£ c t ex is a correlation length for the crumpled-to- 

tubule transition, t y = (T — T ct )/T ct , T ct is the crumpled- 
to-tubule transition temperature, and t y > corresponds 
to the crumpled phase. 

Note that we have built into the scaling laws the fact 
that both R y and Rg scale like L Va in the crumpled 
phase, with v c the radius of gyration exponent for the 
crumpled phase (which, as noted earlier, is the same for 
anisotropic and isotropic membranes). Due to the ex- 
tended nature of the tubule phase, vf = 1, of course. 
The anisotropy in manifested in the crumpled phase only 
through the different temperature dependences of Rg 
and R y . The former of these vanishes as t y — > + (since 
the radius of gyration in the tubule phase is much less 
than that in the crumpled phase, since v% < v c ), which 
implies 7+ > 0, while the latter diverges as t y — ► + , 
since the tubule ultimately extends in that direction, 
which implies 7^ < 0. 



Note also that our expression Eq.7.39, and, in partic- 



ular, the fact that Rg 7^ Ry even above the crumpled- 
to-tubule transition (i.e., in the crumpled phase) , implies 
a spontaneous breaking of rotational invariance even in 
the crumpled phase! This seemingly bizarre (but cor- 
rect) result is actually not all that unfamiliar: polymers, 



as we did earlier in our discussion of mean field theory in 
the absence of self-avoidance, and using p — M/V for the 
embedding space density of the tubule, and again using 
the fact that the tubule mass M ~ L 



D- 



1 L y , we see that 



E SA 



LyL 



2(D-1) 



(y R G 



(7.43) 



Usin g this estimate of the self-avoidance energy in 
Eq.2.1, and estimating the other terms in that expression 
by scaling, we obtain the full Flory theory for the tubule 
phase, with all temperature dependent effects taken (ad- 
mittedly crudely) into account: 



Efl — 



lyCy+UyyCy+tx ( ^ 



LyL 



2(D-1) 



d-1 



Minimizing this over Rg, we obtain 

TV, I U 



Rc 



(y-R-G 
(7.44) 



(7.45) 



where, as we found earlier, v t — ^jj' , but now we have 
the singular temperature dependence of Rg near the 
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crumpled-to-tubule transition explicit through the pres- 
enc e of t he £ y term. Inserting this expression for Rq into 
Eq.[7.44 we find 



EFL = tyCy + Uyy(y + 



», \ 3+T 2{d-D) 

v > - — d+r- \ t d-i , 



L , 



(7.46) 



The exponents defined by Eq.7.39 can now be obtained 
by minimizing Efl in Eq. 7.46| with respect to Q y , which 
amounts to balancing two of the three terms in Efl, 
which two, depending on whether one is interested in the 
crumpled phase (t y > 0), the tubule phase (t y < 0), or 
the transition between them (t y = 0). 

In the crumpled phase t y > 0, as a result, the order 
parameter Q y vanishes in the thermodynamic limit, al- 
lowing us to neglect the quartic C, y term relative to the 
quadratic C, y one. Balancing the remaining two terms 



tyC 



y^y 



, d+l 



(7.47) 



we obtain 



the embedding space. It does so gently by having the 
identical growth (for square membranes) of Rq and R v 
with L, but exhibiting anisotropy via the prefactors, with 
the ratio R v /Rg diverging as the crumpled-to-tubule 
transition is approached. 

The tubule phase is characterized by t y < and a fi- 
nite order parameter £ y > 0. Therefore in this phase, the 
term proportional to t"? m Efl, Eq.7.46, clearly 



becomes negligible relative to the first two terms when 
L± — > oo. Therefore, we can neglect that term for a suffi- 
ciently large membrane (i.e., a membrane larger than the 
critical correlation length £ cr ). Minimizing the remaining 
first two terms in Efl therefore gives C, y oc \/]ty[) (in- 
dependent of Lj_) as in mean-field theory in th e abs ence 
of self- avoidance. Inserting this inside R y , Eq. 7.42| and 
comparing with the general scaling form for R y , implies 
for a square membrane 



«f = 1 



y 

T- = 



(7.54) 
(7.55) 



Using this in the earlier expre ssion Eq. 7.45 for Rq, we 
obtain the last line of Eq. 7.39| , with 



2 ,d-l\2p+2T d _ D 
U -L I T---3+2- 

—d+T L ± 



(7.48) 



Using this expression for £ y inside Eq j7.42| for R y gives 
for a square membrane (L y = L± = L) 



Ru 



1 _ d+l p + 2 

(VH^ 1 ) 2(d+2) ty 2<d + 2) Ll +2 



which, after comparing with the general form for R y , 
Eq.7.39, gives 



D 



d + 2 ' 
d+l 
~2(d + 2) 1 



(7.50) 
(7.51) 



equation 7.50 being a well-known Flory result for the 
radius of gyration exponent v c forj-a. iidimensional man- 
ifold, embedded in d dimensions J23l3 and 7^. new and 
special to anisotropic membranes. Furthermore, insert- 
ing ( y , Eq.7.48 inside Eq. 7.45 for Rq, we obtain 



„G 



D + l 



1 



2(d + l) 



(7.56) 

(7.57) 
(7.58) 



Finally, right at the crumpled-to-tubule transition, 
(7.49) ty = and we must balance the last two terms in Efl, 



Eq. 7.46 . Minimizing Efl over £ y , we find at the transi- 
tion 



Cy OC L 



(d-D) 
3 + 2d 



(7.59) 



which, when inserted in Eq. 7.42 for R y implies for a 
square membrane that 



D + d + 3 

Ry CX L 3 + 2d 



right at the transition. This leads to 
D + d + 3 



,v - 



2d + 3 



(7.60) 



(7.61) 



Rg 



vT d+2W+T)t A _ 2{d+2){d+1) )t 2id+2) Ll +2 , (7.52) 



for a squa re membrane. Using the result Eq. 7.59 for £ y 
in Eq.7.45 for Rq gives, right at the transition, 



which, not surprisingly gives the same expression for v c 
as in Eq.7.50, and predicts 



7? 



1 



2(d+2) ' 



(7.53) 



7+ 7^ 7+ supports our earlier claim that even the crum- 
pled phase spontaneously breaks rotational invariance in 



with 



Rg oc L, 



1 



d-D 
d+l \ J+~2d 
2D + 3 



2d + 3 



(7.62) 

(7.63) 
(7.64) 
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The scaling relations Eqs.7.40 and 7.41, quoted above, 
then give 



2(d — D) 
2d + 3 



(7.65) 



and are reassuringly consistent with our independent cal- 
culatio ns of exponents 7 



G,y 



,G,y 



and z/ ct '* ', Riven 



in Eqs j7i5lL^5^ , ^55|j 7i57L fr750| j7. 54| J7. 5^,[7T6li and f^H 

above. For the physical case of a two dimensional mem- 
brane embedded in a three dimensional space, {D = 
2,d = 3) 



v c = 4/5 , 


(7.66a) 


- 7/9 , 


(7.66b) 


!/» = 8/9 , 


(7.66c) 


"t = 3/4 , 


(7.66d) 


7? - 1/10 , 


(7.66e) 


ll = -2/5 , 


(7.66f) 


7 G = -1/8 , 


(7.66g) 


1 V - = 1/2 • 


(7.66h) 


= 2/9, 


(7.66i) 



Note that the signs of the 7 < f ? / y _ imply that Rq shrinks 
as the crumpled-to-tubule transition is approached from 
above, and grows as it is approached from below, while 
R y does the opposite. Note also that the crumpled-to- 
tubule transition is quite rounded by finite size effects, 
even for large membranes, because of the small value of 
the crossover exponent 4>, which leads to a large correla- 
tion length £ c t{ty)- Taking an example of a L — 10/im 
membrane with lattice constant a = 10A, we find that 
the crumpled-to-tubule transition is rounded at a reduced 
temperature t y ss (L/a)~^ ~ 0.13, while our hypotheti- 
cal simulation of a 10 4 particle net experiences rounding 
at t y « 0.36. Thus, the transition may not appear sharp 
experimentally or in simulations, even though it is, in 
principle, in the thermodynamic limit. 

The singular parts of other thermodynamic vari ables 
obey scaling laws similar to that for Ro,y, Eq. 7.39 . For 
example the singular part of the specific heat per particle 
C v , i.e., a second derivative of the intensive free energy 
with respect to temperature, is given by 



L D dtl \2 v y 



(7.67) 



which, using Eq.7.39 leads to the scaling form for C v 



C v = LPg{t y L*) , 

ty a+ LP- a +*, t v >0,L»£ ct 
L « U 

\ty\- a -LP- a -*, ty < 0,L»£ ct 



(7.68) 



where, 



SO) 



dx 2 y 



(7.69) 



Using the exponents characterizing R y derived above, we 
obtain: 



, Flory theory 



-2 7 * + 1 
2d + 3 



d + 2 



Flory theory 



Flory theory, d = 3 



-2 7 ^ + 1 , 

, Flory theory 



(7.70a) 
(7.70b) 



(7.71a) 
(7.71b) 

(7.71c) 



(7.72a) 
(7.72b) 



This leads to the unusual feature that outside the crit- 
ical regime (i.e. for L » £ ct ), the singular part of 
the specific heat above the crumpled-to-tubule transi- 
tion vanishes in the thermodynamic limit like L~ Q+< ^ ~ 
£-2(d-D)/(d+2) ^ £~ 2 / 5 ; i n the last expression we have 
used the Flory estimates of the exponents, evaluated in 
D = 2 and d = 3. Only within the critical regime does 
the singular part of the specific heat per particle becomes 
nonvanishing as L — > 00. Similar results were first found 
for_t.he direct crumpled-to-flat transition by Paczuski et 
al.Ei 

We now turn to the tubule-to-flat (tf) transition. On 
both sides of this transition, R y = L y x 0(1). Therefore 
only the other two radii of gyration R x and R z exhibit 
critical behavior, which can be summarized by the scaling 
law: 



oc < 



t ± + L 

{ uv- 



t±>0,L » £ t f 
L « fa 
t x <0,L» £ tf 



(7.73) 



where t± = (T — T t f)/T t f, t± > is assumed to corre- 
spond to the tubule phase, £tf oc {bx^ 1 '^** is the corre- 
lation length for this transition, and the exponents obey 
the scaling relations 



uf = C « 0.59 , 
v* f = 1 , 

x.z 

xz V i~ V ti 

7+ - 



7- 



'tf 



(7.74a) 
(7.74b) 

(7.74c) 
(7.74d) 



In the above we have taken the x-direction to be the new 
(in addition to y) extended direction in the flat phase 
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(which is why isf = 1), and £ is the roughness exponential 
of the flat phase (quoted for the physical case D = 2 and 
d = 3), giving the transverse height fluctuations of the 
d — 2 components of the displacement perpendicular to 
the flat membrane. 

To calculate these exponents, we can use Flory the- 
ory in the tubule phase, and at the transition, while in 
the flat phase, where as discussed above, self-avoidanee 
is irrelevant, we simply match onto the scaling theoryeS 
of the flat phase. Doing so, we find that Flory theory 
predicts identical behavior for R x and R z in the tubule 
phase and at the transition: 



D + 3 
d + 3 



= -, for D = 2 , d = 3 , 
6 



i 



i 



d + l ' 

for D = 2 , d = 3 . 



(7.75a) 
(7.75b) 
(7.75c) 
(7.75d) 



We believe that the identical temperature (t±) and scal- 
ing (with L) behavior of R x and R z as the tubule-to-flat 



transition is approached from the tubule side (Eqs. 7.75a 
and 7.75c) is an artifact of Flory theory and that in fact 
R x >> R z throughout this region, with the ratio R x /R z 
actually diverging as the transition is approached from 
above. That is, we expect that in reality vf* > v** and 

7+ < 7+- 

In addition, Flory theory predicts 



4>tf 



2(d-D) 
d + 3 



- , for D = 2 , d = 3 
o 



(7.76) 
(7.77) 



In the flat phase, 7^ follows from simply minimizing 
the mean field free energy without self-avoidance (since 
self- avoidance is irrelevant in the flat phase), giving 



7- 



(7.78) 



while matching R z = -L^i^l 7 - onto the critical pre- 
diction R z cx L Ut f at the correlation length L = £ = 
l*x| 



- 1 /<t>tf gives 



tf 



<t>tf 

-0.73 



(7.79) 
(7.80) 



where the first equality is an exact scaling law, while the 
second, approximate one uses-EJory theory for 4>tf an d 
v^p and the SCSA calculationta of C for the flat phase, 
all evaluated in the physical case D — 2 and d = 3. 

As the tubule-to-flat transition is approached from be- 
low (the flat phase side) R x shrinks as R x « ItjJ 1 ' 2 -^ an d 



R z increases as R z w |i^| -a73 L - 59 with vanishing 
Approaching this transition from above (the tubule phase 
side) R x and R z both extend as R XiZ w |<j_| _1 / 4 L 3 / 4 with 
vanishing tj_ to the L 5 / 6 scaling at the tubule-to-flat crit- 
ical point. 

The singular part of the specific heat again obeys a 
scaling law: 



C v = L^fg tf {txL* t f) , 



cx < 



where, in Flory theory, 



t± > 0,L » fa 
L « fa 
t y < Q,L » £ tf 



(7.81) 



tf 
a_i = 



tf = 2v tf + <f> tf 



(7.82a) 

(7.82b) 
(7.82c) 



Thus, again, the singular part of the specific heat van- 
ishes (now like L -1 / 2 ) in the thermodynamic limit above 
(i.e., on the tubule side) of the transition, while it is 0(1) 
and smooth as a function of temperature in both the crit- 
ical regime and in the flat phase. 



VIII. SUMMARY AND CONCLUSIONS 

In summary, we have studied the effects of intrinsic 
anisotropy in polymerized membrane. We found that this 
seemingly innocuous generalization leads to a wealth of 
new phenomena, most remarkable of which is that any 
amount of anisotropy leads to a new, tubule phase which 
intervenes between the previously predicted flat and 
crumpled phases in anisotropic membranes (See Fig^). 
We have presented a detailed theory of the anisotropic 
membrane focusing on the new tubule phase. Consid- 
ering thermal fluctuations in the tubule phase we have 
shown that the phantom tubule phase exhibits anoma- 
lous elasticity, and calculated the elasticity and size expo- 
nents exactly, as summarized in Eqs. 5. 15, |5. 16 , 5.3 j .5. 
We then considered the physically more relevant case 
of a self-avoiding tubule, finding that self-avoiding in- 



teraction is important for physical dimensionalities. Es- 
tablishing relations between the exponent characterizing 
the diameter of the tubule and the exponents describing 
anomalous elasticity and transverse undulations, we cal- 
culated the tubule diameter, size of the undulations and 
the anomalous elasticity within the Flory and e = d uc —d- 
expansion theories. We have also studied self-avoidance 
within a Gaussian variational approximation, which un- 
fortunately but, we believe incorrectly predicts that self- 
avoiding interaction destroys the tubule phase (as it does 
the crumpled phase) for d < 4. We studied the crumpled- 
to-tubule transition in mean-field theory and with the 
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e = 4 — /^-expansion. Finally we developed a scaling the- 
ory of the crumpled-to-tubule and tubule-to-flat transi- 
tions. 

Our exact predictions for the phantom tubules 
Eqs.5.15 5.16] ,5.38, |5.48 have been quantitatively verified 
in the recent simulations by the authors of Ref. ^. 

The possibility of the existence of a new tubule 
phase intermediate between the fully disordered crum- 
pled phase and fully ordered flat phase is exciting from 
both basic physics and potential applications points of 
view. Recently much attention has focussed on utiliz- 
ing self-assembled microstructures for encapsulations for 
various applications, most notably controlled and slow 
drug delivery.E3 The structural stability of polymerized 
membranes is superior to their liquid membrane analogs. 
The theoretical discovery of the tubule phase significantly 
expands the number of possibilities, and also offers the 
potential tunability (by, e.g., adjusting the strength of 
self-avoidance) of the tubule diameter and therefore the 
amount of encapsulation and rate of delivery. 

The realization of the tubule phase in polymerized 
membranes carries even more significance if the claims 
that the fully crumpled phase in polymerized membranes 
does not exist are in fact correct, since in this case the 
tubule phase is the only disordered phase of a polymer- 
ized membrane. 

With the recent focus on self-assembly, it may be pos- 
sible in the near future to freeze in intrinsic anisotropy by 
polymerizing tilted phase of liquid membranes or cross- 
linking polymers. Further numerical simulations which 
include self-avoidance offer another avenue to investigate 
our predictions. We hope that our work stimulates fur- 
ther theory, simulations and experiments in this area. 



ACKNOWLEDGMENTS 



We both thank The Institute for Theoretical Physics 
at the UCSB, and the organizers of the Biomcmbranes 
Workshop held there, where this work was initiated, for 
our good tans, their hospitality and financial support 
under NSF Grant No. PHY94-07194. We acknowl- 
edge partial support from the Colorado Center for Chaos 
and Complexity through the Summer 1997 Workshop on 
Nucleation and Critical Phenomena in Complex Non- 
linear Systems, during which this work was completed. 
Leo Radzihovsky acknowledges support by the NSF CA- 
REER award, through Grant DMR-9625111, and partial 
support by the A. P. Sloan Foundation. John Toner ac- 
knowledges financial support by the NSF through Grants 
DMR-9625111 and DMR-9634596. 



1 L. Radzihovsky and J. Toner, Phys. Rev. Lett. 75, 4752 
(1995). 

2 M. Bowick, M. Falcioni, and G. Thorleifsson, Phys. Rev. 



Lett. 79, 885 (1997); |cond-mat/9705059| . 
3 M. Bowick an d E. Guitter, Phys. Rev. £56, 1 (1997); |cond- 



mat/9705045 



For a review, see the articles in Statistical Mechanics of 
Membranes and Interfaces, edited by D. R. Nelson, T. Pi- 
ran, and S. Weinberg (World Scientific, Singapore, 1989). 

5 D. R. Nelson and L. Peliti, J. Phys. (Paris) 48, 1085 (1987). 

6 P. Hohenberg, Phys. Rev. 158, 383 (1967); N. D. Mermin 
and H. Wagner, Phys. Rev. Lett. 17, 1133 (1966); S. Cole- 
man, Commun. Math. Phys. 31, 259 (1973); 

7 J. A. Aronovitz and T. C. Lubensky, Phys. Rev. Lett. 
60, 2634 (1988); J. A. Aronovitz, L Golubovic, and 
T. C. Lubensky, J. Phys. (France) 50, 609 (1989). 

8 E. Guitter, F. David, S. Leibler, and L. Peliti, Phys. Rev. 
Lett. 61, 2949 (1988). 

9 L. Radzihovsky, unpublished. 

10 A. Rudolph, J. Calvert, P. Schoen, and J. Schnur, "Tech- 
nological Developments of Lipid Based Tubule Microstruc- 
tures" , in Biotechnological Applications of lipid microstruc- 
tures, ed. B. Gaber et al. (Plenum, 1988); R. Lipkin, Sci- 
ence, 246, 44 (December 1989). 

11 Actually, in most experimental realizations, the polymer- 
ization is random and leads to an isotropic disordered mem- 
brane, making the isotropic case quite general. 

12 E. Sachmann, private communication. 

13 D. Bensimon, private communication. 

14 M. Doi and S. F. Edwards, The Theory of Polymer Dynam- 
ics, Oxford University, New York, 1986; P. G. de Gennes, 
Scaling Concepts in Polymer Physics, Cornell University, 
Ithaca, 1979. 

15 G. Grinstein and R. Pelcovits, Phys. Rev. Lett AT , 856 
(1981); Phys. Rev. AA26, 915 (1982). 

16 In the case of polymers, Flory theory agrees with the ex- 
act predictions for the radius of gyration exponent v in 
all dimensions d where such exact predictions exist; in 
d = 4,d = 2, and d = 1, Flory theory recovers the exact 
results of v = 1/2, v — 3/4, and v = 1, respectively. And 
in d = 3 dimensions (where an exact result is not available) 
it agrees with the e-expansion to better than 1%. 

17 G. Grest, J. Phys. I (France), 1, 1695 (1991). 

18 M. Plischke and D. Boal, Phys. Rev. A 38, 4943 (1988); 
F. F. Abraham, W. E. Rudge, and M. Plischke, Phys. Rev. 
Lett 62, 1757 (1989); J. S. Ho and A. Baumgartner, Phys. 
Rev. Lett. 63, 1324 (1989); F. F. Abraham and D. R. Nel- 
son, J. Phys. France 51 2653 (1990); Science 249, 393 

(1990) . 

19 A. Baumgartner and W. Renz, Europhys. Lett. 17, 381 

(1991) ; A. Baumgartner, J. Phys. I (France) 1, 1549 
(1991). 

20 The failure to find the crumpled phase in computer simu- 
lations for self-avoiding membranes might be because self- 
avoidance interaction shifts the bare k to high values, such 
that in the almost universally employed numerical ball and 
spring models this always put the isotropic membrane on 
the flat side of the crumpling transition. 

21 M. Goulian, J. Phys. II (France), 1, 1327 (1991). 

22 P. Le Doussal, J. Phys. A, 25, L469 (1992). 



33 



23 This is analogous to the familiar 1/n expansion for criti- 
cal phenomena, in which one expands about the number of 
spin components n — » oo limit. 

24 M. Paczuski, M. Kardar and D.R. Nelson, Phys. Rev. Lett. 
60, 2638 (1988). 

25 M. E. Fisher and D. R. Nelson, Phys. Rev. Lett, 32, 1350 
(1974); D. R. Nelson, J. M. Kosterlitz, and M. E. Fisher, 
Phys. Rev. Lett, 33, 813 (1974), and Phys. Rev. B13, 
412 (1976); A. Aharony, in Phase Transitions and Critical 
Phenomena , edited by C. Domb and M. S. Green (Aca- 
demic, New York, 1976), Vol. 6, and Bull. Am. Phys. Soc. 
20, 16 (1975). 

26 J. Toner, Phys. Rev. Lett. 62, 905 (1988). 

27 L.D. Landau and E.M. Lifshitz, Theory of Elasticity, Perg- 
amon Press, (1975). 

28 P. Le Doussal and L. Radzihovsky, Phys. Rev. Lett. 69, 
1209 (1992). SCSA is incredibly successful in that for the 
flat phase of polymerized membranes it predicts exponents 
that are exact in d — > oo, d — D and correct to a leading 
order in e = 4 — D, thereby showing agreement with all 
known exact results. 

29 M. Falcioni, M. Bowick, E. Guitter, and G. Thorleifsson, 
Europhys. Lett. 38, 67 (1997). 

30 M. Kardar and D. R. Nelson, Phys. Rev. Lett. 58, 1289 
(1987); 58, 2280(E) -Phys. Rev. A 38, 966 (1988). 

31 Y. Kantor and D. R. Nelson, Phys. Rev. Lett. 58, 
1289 (1987); Y. Kantor, M. Kardar, and D. R. Nelson, 
Phys. Rev. A 35, 3056 (1987). 

32 J. A. Aronowitz and T. C. Lubensky, Europhys. Lett. 4, 395 



(1987). 

' L. Radzihovsky and J. Toner, unpublished. 
In doing this power counting, o ne m ust assume that the 
q = value of the integral in Eq. 5.13 is precisely cancelled 
by the bare g y . This assumption is the only way to solve 
Eq|OJ for D < D uc . 

B. Nienhuis, Phys. Rev. Lett. 49, 1062 (1982); B. Du- 
plantier in Fields, Strings and Critical Phenomena, Les 
Houches Lectures, edited by E. Brezin and J. Zinn- Justin 
(1990). 

' K. G. Wilson and J. Kogut, Phys. Rep. C, 12, 77 (1977). 
S. F. Edwards, unpublished 1968; J. des Cloizeaux, Kyoto 
Conference on Statistical Mechanics, J. Phys. Soc. Japan 
Suppl. 26, 42 (1969); J. Phys. France 31, 715 (1970). J. des 
Cloizeaux and G. Jannink, Polymers in Solution: Their 
Modelling and Structure (Oxford University Press, Oxford, 
1989). 

' We thank Stephanie Palmer who generously donated her 
time to us (computer illiterates), t o num erically evaluate 



7.32 



and 7.37 



integrals appearing in Eqs. 

' R. P. Feynman, Statistical Mechanics: A Set of Lectures, 
Addison- Wesley (1972). 

1 T. A. Tokuyasu and J. Toner, Phys. Rev. Lett. 68, 3721 
(1992). 

B. Jacobsen and L. Radzihovsky, work in progress. 
L. Cannon, J. Aronovitz, and P. Goldbart, J. Physique J, 
1, 629 (1991); J. Aronovitz and M. Stephens, J. Phys. A 
20, 2539 (1987); J. Aronovitz and D. Nelson, J. Physique 
47, 1445 (1986). 



34 



